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1 .  INTRODUCTION 


1 . 1  Summary 

Continuous  fiber  composite  laminates  are  known  to  undergo  a 
substantial  amount  of  complex  load-induced  damage  which  can  adversely 
affect  component  performance  C 1 3 -  Therefore,  it  is  desirable  to  develop 
new  models  capable  of  accounting  for  the  effect  of  damage  on  materials 
properties. 

This  report  documents  research  completed  during  the  first  year  of 
a  three  year  effort  under  AFOSR  grant  no.  AFOSR-84-0067  and  originally 
detailed  under  Texas  A&M  Research  Foundation  proposal  no.  RF-84-3M  and 


dated  October  1983.  The  objective  of  this  research  is 

to  develop  an 

accurate 

damage  model  for  predicting  strength  and 

stiffness  of 

continuous 

fiber  composite  media  subjected  to  fatigue 

or  monotonic 

loading  and  to  verify  this  model  with  experimental  results 

obtained  from 

composite 

specimens  of  selected  geometry  and  makeup  to 

be  described 

herein. 


1  . 2  Statement  of  Work 

The  following  is  a  brief  summary  of  work  to  be  performed  under  the 
present  grant: 

1)  develop  constitutive  equations  relating  stresses  to  strains  and 
damage  internal  state  variables  (ISV)  which  may  be  used  in  a  stress 
gradient  field; 

2)  develop  ISV  growth  laws  as  a  function  of  load  history  for 
matrix  cracking,  interlaminar  fracture,  etc.; 


3)  develop  finite  element  algorithms  capable  of  evaluating  ply 
properties  in  damaged  components. 

*1)  perform  experiments  on  components  with  selected  stacking 
sequences  in  order  to  verify  the  model. 

2.  RESEARCH  COMPLETED  TO  DATE 

2. 1  Summary  of  Completed  Research 

The  following  research  has  been  completed  during  the  first  year: 

1)  the  general  damage-dependent  stress-strain  relations  have  been 
revised  and  completed; 

2)  an  ISV  growth  law  has  been  constructed  for  matrix  cracking; 

3)  laminate  equations  have  been  constructed  for  matrix  cracking; 
the  relation  between  the  damage  ISV  and  surface  area  of  cracks 

has  been  established  analytically; 

5)  all  finite  element  programs  are  complete;  and 

6)  the  model  has  been  compared  to  both  experimental  and  finite 
element  results  for  glass/epoxy  laminates. 

In  addition,  the  following  research  is  partially  completed; 

1 )  experimental  characterization  of  ply  and  laminate  properties  in 
a  variety  of  graphite/epoxy  ASil/3502  laminates; 

2)  non-destructive  evaluation  of  damage  in  various  graphite/epoxy 
laminates;  and 

3)  correlation  of  model  to  experiment  for  degraded  stiffness  due 
to  damage  in  graphite/epoxy  composites. 

The  following  sections  further  detail  the  results  summarized 


above. 
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2.2  Development  of  the  General  Model 

Thermodynamic  and  symmetry  constraints  have  been  utilized  to 
construct  stress-strain  relations  of  the  form: 
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where  a  is  the  stress  tensor,  e  is  the  strain  tensor,  o' 
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damage  dependent  residual  stress,  e.  is  the  damage  dependent  thermal 
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strain,  and  C'  is  the  damage  dependent  modulus  tensor,  given  by 
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where  C  and  M.n..  ,  are  material  constants  and  a,n  are  second-order 
L .  . ,  ,  ljklmn  u 

ijkl  mn 

tensor- valued  damage  dependent  internal  state  variables. 

A  kinematically  based  description  of  the  internal  state  variable  has 


been  constructed  of  tha  form 
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where  u.  are  the  components  of  the  crack-opening  displacement,  n^  are 

the  components  of  the  crack  unit  normal,  and  is  the  total  surface 

area  of  cracks  in  a  given  volume  element,  V^,  as  shown  in  Fig.  1. 

An  ISV  growth  law  has  been  constructed  for  matrix  cracking,  as 

given  by 
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where  the  parameters  are  as  described  in  Appendix  6.2. 


The  results  outlined  in  this  section  are  described  in  detail  in 


Appendix  6.1. 


2.3  Application  to  Matrix  Cracking 

The  model  has  been  applied  to  composite  laminates  with  a  single 
damage  mode  consisting  of  matrix  cracks,  as  shown  in  Fig.  2.  This  was 
accomplished  by  imposing  symmetry  constraints,  reducing  to  single  index 
notation,  and  constructing  laminate  equations.  It  was  found  that  for 
matrix  cracks  which  lay  in  a  plane  normal  to  the  fiber  (x^  direction, 
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in  local  ply  coordinates,  where  a 


is  the  damage  tensor  for  a  typical 


ply  due  to  matrix  cracking. 

An  examination  of  the  dependence  of  free  energy,  u^,  on  damage  led 

to  the  conclusion  that  the  magnitude  of  damage  |  |  depends  only  on  the 

surface  area  of  cracks  S,,.  Furthermore,  since  the  free  energy  depends 

l  1| 

explicitly  on  damage,  it  follows  that  the  relation  between  la^l  and  S2 
could  be  determined  from 
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where  G  is  the  energy  release  rate. 
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(6) 


The  energy  release  rate  was  then  calculated  as  a  function  of  by 
utilizing  a  finite  element  simulation  of  a  [0,90]s  laminate. 
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Crack  Face 
Tangent 


Fig.  2.  Description  of  the  Internal  State  at  Point  O'. 


Using  the  procedure  briefly  described  above,  it  is  possible  to 


relate  all  stiffness  components  to  cracked  surface  area  for  composite 
laminates  with  matrix  cracks.  It  is  emphasized  that  the  model  predicts 
all  stiffness  components  for  any  layup  using  material  property  input 
which  is  stacking  sequence  independent. 

Figs.  3  &  4  show  predictions  of  the  model  for  axial  and  other 
stiffness  components  of  [0,90-]  E-glass  epoxy  laminates.  Fig.  3 

J  5 

indicates  that  the  model  predicts  somewhat  higher  stiffness  than 

experimental  results.  Although  experimental  results  are  not  available 

for  other  components,  Fig.  M  shows  excellent  agreement  between  FEM 

results  and  the  model  for  v  .  Details  of  the  results  outlined  in  this 

xz 

section  are  contained  in  Appendix  6.2. 

The  following  section  will  detail  attempts  to  obtain  more  useful 
experimental  results. 


2.M  Experimental  Research  Activities 

I  The  purpose  of  the  experimental  research  is  to  study  and  document 

the  progression  of  damage  in  a  systematic  manner  that  will  fully  support 
the  development  of  the  constitutive  model.  The  authors  have  determined 
via  an  extensive  literature  survey,  partially  documented  in  reference  1 , 
that  the  experimental  data  necessary  to  characterize  and  verify  the 
model  does  not  exist  in  the  current  literature.  The  comprehensive 
experimental  data  base  will  provide  two  fundamental  types  of  essential 
information  which  are  not  currently  available  in  the  open  literature. 
First,  the  progression-of-damage  study  will  establish  the 
phenomenological  characterization  of  the  internal  state  variables  and 
will  provide  the  essential  ingredients  for  formulating  the  associated 
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damage  growth  laws.  Second,  the  experimental  data  base  will  provide 
direct  measurements  of  the  effective  (damage-degraded)  moduli.  These 
direct  measurements  will  be  compared  to  the  constitutive  model 
predictions.  This  is  an  essential  task  for  both  the  development  and 
"fine-tuning"  of  the  model  as  well  as  establishing  the  validity  and 
accuracy  of  the  model.  The  following  sections  describe  the  results  of 
the  three  primary  tasks  that  have  been  completed  to  date.  The  first 
task  was  the  development  of  a  materials  character ization  laboratory. 
The  second  task  was  the  formulation  of  a  test  matrix,  laminate  panel 
preparation  and  specimen  coupon  fabrication.  Finally,  the  study  of 
matrix  crack  damage  has  commenced. 

2.4.1  Materials  Characterization  Laboratory 

The  research  objectives  of  this  grant  necessitated  the 
establishement  of  a  materials  characterization  laboratory.  The  effort 
was  supported  jointly  by  Texas  A&M  University  and  by  an  equipment  grant 
from  AFOSR.  The  established  laboratory  meets  the  dual  requirements  of 
measuring  the  mechanical  properties  of  materials  and  nondestructi vely 
studying  load-induced  microstructural  damage.  The  primary  laboratory 
facilities  consist  of  a  20-kip  Instron  testing  machine  for  monotonic 
loading  conditions  and  a  110-kip  MTS  880  testing  system  (partially 
funded  by  AFOSR  Equipment  Grant-84-0257)  for  both  monotonic  and  cyclic 
loading  conditions.  Support  equipment  include  computer  controlled  data 
acquisition  and  data  reduction  systems,  as  well  as  redundant  real-time 
analog  signal  recording  systems.  (The  MTS  880  testing  system  also 
includes  an  environmental  chamber  for  cryogenic  or  elevated  temperature 
testing  in  a  controlled  humidity  environment.)  The  primary 
nondestructive  evaluation  (NDE)  facility  is  a  portable  x-ray  radiography 


system  for  in-situ  x-ray  examinations.  A  second  NDE  facility  is  an 
ultrasonic  spectroscopy  laboratory.  This  technology  is  being  developed 
by  Dr.  V.K.  Kinra  as  part  of  a  "sister"  AFOSR  Grant.  Additional  NDE 
facilities,  already  in  place  at  Texas  A&M  include  edge-replication, 
ultrasonic  C-scan,  optical  microscopes  and  a  scanning  electron 
microscope.  All  of  the  above  equipment  is  now  in  place  and  functional 
along  with  fully  trained  personnel. 

2.4.2  EXPERIMENTAL  TEST  PROGRAM 

Laminate  panels  are  being  prepared  from  pre-impregnated 

l 

graphite/epoxy  tape,  AS4/3502,  by  the  Mechanics  and  Materials  Laboratory 
at  Texas  A&M  University.  The  curing  process  was  developed  from  the 
procedure  recommended  by  the  pre-preg  tape  vendor.  Each  12"  x  12"  panel 
yields  ten  1"  x  11"  tensile  test  coupons.  The  matrix  of  test  coupons  is 
given  in  Table  1.  The  anticipated  specimen  requirements  necessitated 
the  fabrication  of  2  panels  for  each  laminate  stacking  sequence.  While 
panel  fabrication  is  still  an  ongoing  activity,  all  Type  I  and  II 
panels  have  been  fabricated.  The  cured  materials  are  being  stored  in  a 
dessicant  chamber. 

The  number  of  specimens  in  each  test  category  listed  in  Table  1 
were  selected  primarily  to  insure  adequate  replicate  test  data.  The 
specified  load  histories  for  the  four  fatigue  tests  are  considered  to  be 
tentative.  The  exact  cyclic  load  histories  will  be  selected  at  a  later 
date.  For  this  reason,  it  is  desirable  to  have  a  sufficient  number  of 
specimens  available  should  the  research  warrant  a  redirection.  It 
should  be  noted  that  there  are  two  spare  specimens  in  addition  to  the 
eighteen  specimens  listed  in  the  test  matrix  for  each  laminate.  The 


TABLE  1  TEST  MATRIX 
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B.  LAMINATE  CHARACTERIZATION  &  DAMAGE  GROWTH  LAWS 
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total  number  of  specimens  (20  per  laminate)  will  allow  for  some 
destructive  examinations  of  the  microstructural  damage  without 
diminishing  the  overall  data  base. 

The  laminates  listed  in  Table  1  were  selected  for  specific 
reasons.  The  Type  I  laminates  were  selected  to  study  matrix  cracks. 
The  primary  mode  of  microstructural  damage  in  cross-ply  laminates  is 
matrix  cracks  in  the  90°  plies  [2,3]  In  addition,  the  free  edge 
interlaminar  stresses  are  minimal  and  should  not  affect  the 
microstructural  damage  that  develops  prior  to  fracture.  Therefore,  the 
cross-ply  laminates  are  ideal  for  isolating  the  single  damage  mode  of 
matrix  cracking.  The  specific  cross-ply  stacking  sequences  were 

selected  so  that  both  the  thicknesses  of  the  90°  ply  layers  and  the  0° 
ply  constraint  layers  are  varied. 

The  quasi-isotropic  laminates,  Type  II,  were  selected  because  of 
their  practical  significance  and  in  order  to  study  adjacent  ply  matrix 
crack  interaction.  The  microstructural  damage  should  be  free  of  edge 
effects  in  both  stacking  sequences  because  they  both  have  compressive 
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interlaminar  normal  stresses  at  the  free  edge  [^].  Also,  the  +  45°  and 
-45°  plies  are  adjacent  in  one  laminate  and  separated  in  the  other. 
This  will  provide  information  concerning  the  effect  of  adjacent  ply 
constraint  on  matrix  cracking  and  crack  interaction.  Finally,  these  two 
laminates  will  be  useful  for  the  study  of  localized  internal 
delaminations  which  accompany  the  intersection  of  adjacent  ply  matrix 
cracks . 

The  Type  III  laminates  were  chosen  to  isolate  and  study  matrix 
crack  damage  in  (+)  and  (-)  ;■  plies  in  the  absence  of  matrix  cracks  in 
90°  plies.  The  two  laminate  stacking  sequences  are  consistent  with  the 
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Type  II  laminates,  excluding  the  90°  plies,  and  both  have  minimal 

interlaminar  normal  stress  at  the  straight  free  edges. 

The  Type  IV  laminates  were  selected  to  study  the  influence  of  free 

edge  effects  [2,5].  The  failure  of  the  [ 302/— 302  3S  laminate  is 

dominated  by  the  interlaminar  shear  stress,  t  .  Free  edge  delaminations 

xz 

open  at  the  +30/-30  interface.  On  the  other  hand,  delaminations  do  not 
form  at  the  free  edges  of  the  [(±30)2^  laminate.  The  significance  of 
free  edge  effects  is  illustrated  by  the  fact  that  the  strength  of  the 
[(±30)2^  laminate  is  much  higher  than  that  of  the  L302<'~30^3  laminate. 
It  is  anticipated  that  free  edge  delaminations  can  only  be  addressed  as 
a  boundary  value  problem  since  the  delamination  surfaces  form  external 
boundaries.  However,  this  damage  mode  is  included  herein  because  it 
represents  a  realistic  damage  mode  for  many  practical  structural 
geometries. 

2.4.3  Experimental  Results 

The  first  phase  of  the  research  has  the  objective  of  establishing 
the  internal  state  variable  and  the  associated  damage  growth  law  for  the 
isolated  damage  mode  of  matrix  cracks.  Type  I  cross-ply  laminates  have 
been  selected  for  this  purpose.  This  initial  phase  of  experimental 
research  is  being  confined  to  the  quasi-static,  monotonic  loading 
condition.  The  general  expeimental  procedure  is  briefly  outlined  below: 

1.  Establish  the  basic  lamina  properties  for  input  to  all  analyses. 

2.  Determine  the  undamaged  moduli  and  the  load-to-failure  versus 
displacement  characteristics  of  the  laminate. 

3.  Using  the  results  of  2  as  a  guide,  study  and  document  the 
progression  of  damage  as  follows: 
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(a)  At  appropriate  load  increments,  perform  nondestructive 
examinations. 

(b)  Determine  the  tangent  moduli  and  secant  moduli.  Record  all 
extensometer  versus  load  data. 

(c)  At  a  few  specially  selected  damage  states,  perform  destructive 
examinations . 

The  basic  lamina  properties  are  listed  in  Table  2.  In  each  case, 
the  mean  value  and  standard  deviation  are  listed.  From  the  limited 
comparisons  available,  it  can  be  seen  that  the  recently  measured  lamina 
properties  are  in  good  agreement  with  the  specified  values  supplied  by 
the  pre-preg  tape  vendor. 

The  experimental  investigation  of  matrix  cracks  in  the  90°  plies 

of  cross-ply  laminates  has  commenced.  A  typical  progr ess ion-of- damage 

pattern  li  the  [0/90.J  laminate  is  shown  in  Fig.  5.  Also,  Fig.  6  shows 

j  s 

an  associated  x-ray  radiograph  of  the  apparent  saturation  damage  state 

in  the  [0/90  ]  laminate.  The  edge  replicas  shown  in  Fig.  5  were  taken 
d  s 

of  the  same  specimen  at  increasing  applied  stress  levels.  The 
increasing  number  of  cracks  and  the  apparent  saturation  spacing  at  67.1 
ksi  are  consistent  with  well  known  results  for  transverse  matrix 
cracking.  It  is  obvious  from  Fig.  5  that  all  of  the  matrix  cracks  are 
not  straight  cracks  through  the  thickness  of  the  90°  layer.  The  initial 
cracks  that  form  are  straight.  However,  as  cracks  begins  to  fill  in  the 
spaces  between  the  initial  cracks,  some  cracks  have  distinctive 
parabolic  profiles.  This  can  clearly  be  seen  in  the  edge  replica  taken 
at  56.8  ksi.  The  curved  cracks  have  greater  surface  area  than  straight 
cracks  and  unit  normals  that  are  not  aligned  with  the  x^axis  direction 
as  is  the  case  with  straight  cracks.  Therefore,  the  presence  of  the 


TABLE  2  MEASURED  LAMINA  PROPERTIES 

MEAN  STANDARD  VENDOR 

VALUE  DEVIATION  SUPPLIED 

PROPERTIES 

21.5xl06  psi  2.0%  21 . 5xl06  psi 

1.39xl06  psi  2.1% 

0.694xl06  psi 

l 

0.31  3.7% 

326.0  ksi  3.5%  310.0  ksi 

0.0144  4.6 % 

1  f  tensile  strength  of  0°  unidirectional  laminate, 

tu^ 

O 

2.  e.  =  tensile  failure  strain  of  0  unidirectional  laminate 
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curved  cracks  may  have  a  significant  influence  on  the  formulation  of  the 
internal  state  variable  for  matrix  cracking.  The  initial  indication  is 
that  the  profile  of  the  curved  cracks  is  greatly  influenced  by  the 
thickness  of  the  90°  layer  and  to  some  extent  by  the  thickness  of  the 
adjacent  0°  plies.  This  is  probably  related  to  the  previously 
identified  shear  lag  effect  on  matrix  crack  formation.  This  is 
currently  the  focus  of  the  experimental  research  effort. 


2.5  Conclusion 


A  general  model  has  been  developed  for  predicting  the  relation 

between  stresses  and  strains  in  composites  with  load-induced  damage. 

This  model  is  capable  of  predicting  structural  component  response  in  a 

stress  gradient  field.  Although  the  general  form  of  the  model  is  fairly 

complex,  it  has  been  shown  herein  that  for  the  case  of  matrix  cracking 

the  model  is  only  slightly  more  cumbersome  that  standard  laminate 

analysis.  Furthermore,  the  model  has  been  compared  with  some  success  to 

experimental  results  for  [0,90  ]  E-glass  epoxy  laminates. 

j  a 

Current  efforts  deal  with  verification  of  the  model  for  ASiJ/3502 
graphite-epoxy  laminates,  with  primary  emphasis  on  matrix  cracking. 
Since  necessary  experimental  results  are  not  available  in  the 
literature,  a  comprehensive  experimental  program  has  been  initiated 
under  the  current  grant.  Initial  results  indicate  that  the  model  is 
accurate  for  graphite-epoxy  systems. 

Activities  in  the  second  year  of  research  will  concentrate  on 
three  important  issues:  1)  effect  of  matrix  cracking  on  reduction  of  all 
stiffness  components  as  a  function  of  stacking  sequence  in  cross-ply 
laminates;  2)  refinement  of  the  ISV  growth  law  for  matrix  cracking;  and 
3)  initial  studies  of  the  effects  of  interlaminar  delamination  on 
stiffness  loss. 

The  authors  are  quite  encouraged  by  the  current  results  and 
believe  that  they  provide  ample  justification  for  continuing  the 


research  effort. 


1  VrsT 


-r,  v-.'  iv'-'.".'*; . 


-T^T 


ri  U  iiin.n 


2 . 6  References 

[1]  Groves,  S.E.,  and  Allen,  D.H.  ,  "A  Survey  of  Damage  in  Continuous 
Fiber  Composites,"  Texas  A&M  University  Mechanics  and  Materials  Center, 
MM-5023-84-6,  March,  1984. 

[2]  Herakovich,  C.T.,  Nagarkar,  A.,  and  O'Brien,  D.A.,  "Failure  Analysis 
of  Composite  Laminates  with  Free  Edges,"  Modern  Developments  in  Comp. 
Mat,  and  Struc.  Edited  by  J.R.  Vinson,  ASME,  345  East  47th,  St.,  NY,  NY, 
10017,  1979,  pp.  53-66. 

[3]  Jones,  R.M.,  Mechanics  of  Composite  Materials,  McGraw-Hill  Book  Co., 
NY,  NY,  1975. 

[4]  Herakovich,  C.T.,  "On  the  Relationship  Between  Engineering 
Properties  and  Delamination  of  Composite  Materials,"  J.  of  Comp.  Mat., 
Vol .  15,  July  1981,  pp.  336-348. 

[5]  Herakovich,  C.T.,  "Influence  of  Layer  Thickness  on  the  Strength  of 
Angle-Ply  Laminates,"  J.  of  Comp.  Mat . ,  Vol.  16,  May  1982,  pp.  216-227. 


3.  PUBLICATIONS  LIST 


3.1.  To  Be  Submitted  Tor  Publication 

1.  "A  Damage  Model  for  Continuous  Fiber  Composites,  Part  I- 
Theoretioal  Development,"  by  D.H.  Allen,  S.E.  Groves,  and  R.A.  Schapery 
(Appendix  6.1),  to  be  submitted  to  Mechanics  of  Materials. 

2.  "A  Damage  Model  for  Continuous  Fiber  Composites,  Part  II  - 
Model  Applications,"  by  S.E.  Groves,  D.H.  Allen,  C.E.  Harris,  R.A. 
Schapery,  (Appendix  6.2),  to  be  submitted  to  Mechanics  of  Materials. 

i 

3.2.  In  Preparation 


1.  "Modelling  of  Stiffness  Reduction  Due  to  Matrix  Cracks  in 
Graphite/Epoxy  Laminates,"  by  C.E.  Harris,  R.G.  Norvell,  D.H.  Allen,  and 
S.E.  Groves. 

2.  "Experimental/Analytical  Correlation  of  Damage  Accumulation  in 
Laminated  Composites,"  R.G.  Norvell,  C.E.  Harris,  D.H.  Allen,  and  I. 
Georgiou. 

3.  "Damage  Growth  Laws  for  Matrix  Cracks  in  Graphite/Epoxy 
Laminates,"  by  D.H.  Allen,  S.E.  Groves,  I.  Georgiou,  and  C.E.  Harris. 


4.2  Additional  Staff 


1.  Mr.  S.E.  Groves  (Lecturer,  Research  Assistant,  and  Ph.D.  Candidate)- 
ISV  growth  laws;  laminate  analysis;  finite  element  modelling;  mechanical 
testing. 

2.  Mr.  R.G.  Norvell  (Research  Assistant  and  M.S.  Candidate)-  mechanical 
testing;  nondestructive  evaluation. 

3.  Mr.  I.  Georgiou  (Research  Assistant  and  M.S.  Candidate)-  mechanical 
testing;  nondestructive  evaluation;  ISV  growth  laws. 

4.  Mr.  B.  Harbert  (Lab  Technician)-  experimental  lab  support. 

5.  Mr.  C.  Fredericksen  (Lab  Technician)-  experimental  lab  support. 

6.  Mrs.  T.  Marquez  (Secretary)-  secretarial  support. 


5.  INTERACTIONS 


5. 1  Papers  Presented 


5.1.1  Presented  Under  This  Grant 

1.  D.H.  Allen,  "A  Damage  Model  for  Continuous  Fiber  Composites," 
Tenth  Annual  Composites  Review,  Dayton,  October,  1984. 

2.  D.H.  Allen,  "A  Damage  Model  for  Continuous  Fiber  Composites," 
Society  of  Engineering  Science  21st  Annual  Meeting,  Blacksburg,  VA, 

l 

October,  1984. 

5.1.2  Related  Research  Presentations 

1.  D.H.  Allen,  "On  the  Use  of  Internal  State  Variables  in 
Thermoviscoplastic  Constitutive  Equations,"  2nd  Symposium  on  Nonlinear 
Constitutive  Relations  for  High  Temperature  Applications,  Clevelend, 
June,  1984. 

2.  W.E.  Haisler,  "Numerical  Considerations  in  the  Development  and 
Implementation  of  Constitutive  Models,"  2nd  Symposium  on  Nonlinear 
Constitutive  Relations  for  High  Temperature  Applications,  Cleveland, 
June,  1984. 

3.  C.E.  Harris,  "Comparison  of  the  Fracture  Behavior  of  Thick  Laminated 
Composites  Utilizing  the  Compact  Tension,  Three-Point  Bend  and  Center- 
Cracked  Tension  Specimens,"  presented  at  the  Seventeenth  Symposium  on 
Fracture  Mechanics  sponsored  by  ASTM,  Albany,  NY,  August  1984. 


4.  C.E.  Harris,  "A  Damage  Tolerant  Design  Parameter  for  Graphite/Epoxy 
Laminated  Composites,"  presented  at  the  21st  Annual  Meeting  of  the 
Society  of  Engineering  Science,  Blacksburg,  VA,  October  15-17,  1984. 

5.  C.E.  Harris,  "The  Fracture  of  Thick  Graphite/Epoxy  Laminates  with 
Part-Through  Surface  Flaws,"  presented  at  the  Symposium  on  Composite 
Materials:  Fatigue  and  Fracture,  Dallas,  TX,  October  1984. 

6.  C.E.  Harris,  "Unique  Aspects  of  Conducting  Strength  Tests  of  Thick 
Composites"  invited  paper  presented  at  the  fall  meeting  of  the  Society 
for  Experimental  Mechanics,  Milwaukee,  WI . ,  November  1984., 

5.1.3  To  Be  Presented 

1.  D.H.  Allen,  "Modelling  of  Stiffness  Reduction  Due  to  Matrix  Cracks 
in  Graphite/Epoxy  Laminates,"  Society  of  Engineering  Science  22nd  Annual 
Meeting,  State  College,  PA,  October,  1985. 

2.  C.E.  Harris,  "Experimental/Analytical  Correlation  of  Damage 
Accumulation  in  Laminated  Composites,"  Society  of  Engineering  Science 
22nd  Annual  Meeting,  State  College,  PA,  October,  1985. 

5.2  Research  Related  Travel  and  Consultative  Functions 

1.  D.H.  Allen  visited  Hercules  Corp.,  Magna,  UT,  to  discuss  composite 


materials,  April  1984. 


2.  D.H.  Allen  visited  Martin-Marietta,  Denver,  to  discuss  composites  in 
large  space  structures  and  to  see  the  PACOSS  experiment,  April  1984. 

3.  D.H.  Allen  and  W.E.  Haisler  attended  the  2nd  AFOSR  Forum  on  Large 
Space  Structures,  Washington,  D.C.,  June  1984. 

4.  D.H.  Allen  and  W.E.  Haisler  attended  the  2nd  Symposium  on  Nonlinear 
Constitutive  Equations  for  Elevated  Temperature  Applications,  Cleveland, 
June  1984. 

I 

5.  D.H.  Allen  attended  NASA  Shuttle  Launch  41 -C  to  see  the  launch  of 
LDEF,  which  contains  one  of  Dr.  Allen's  experiments,  April  1984. 

6.  W.E.  Haisler  attended  the  AFFDL  Vibration  Damping  Conference,  Long 
Beach,  CA,  April  1984. 

7.  W.E.  Haisler  attended  the  AIAA  SDM  Conference,  Palm  Springs,  CA,  May 
1  984. 

8.  D.H.  Allen  consulted  with  engineers  at  the  General  Electric  Aircraft 
Engine  Business  Group  on  viscoplasticity  in  metals  at  elevated 
temperature,  Cincinnati,  April  and  November,  1984. 

9.  D.H.  Allen  visited  Drs.  G.  Sendeckyj  and  D.  Paul  at  AFML  to  discuss 
composites,  Dayton,  November,  1984. 

2  6 


Mechanics  and  Materials  Center 
TEXAS  A&M  UNIVERSITY 
College  Station,  Texas 


A  DAMAGE  MODEL  FOR  CONTINUOUS  FIBER  COMPOSITES 
PART  I:  Theoretical  Development 


..  i 

»  « 


■  i 


D.  H.  Allen 


S.  E.  Groves 


Aerospace  Engineering  Department 


R.  A.  Schapery 

Civil  and  Aerospace  Engineering  Departments 
Texas  A&M  University 
College  Station,  Texas  77843 


MM-5023-84-17 


August  1984 
Revised  February  1985 


I»4«» 


REPORT  DOCUMENTATION  PAGE 


la  REPORT  SI  L^HiTv  C  LASSlf  LAI  ION 

Unclassif ied 

lb  RESTRICTIVE  MARKINGS 

2a  SECURITY  CLASSIFICATION  AuTHOHiTY 

3  D'STRIBUTION/AVAILABILITY  of  REPORT 

2b.  OECLASSIF  iCAT  ION,  DOWNGRADING  SCHEDULE 

Unlimited 

4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

_ _ 1 

5.  monitoring  organization  report  numberisi 

6a.  NAME  CF  PERFORMING  ORGANIZATION 

Aerospace  Engineering  Dept. 

6b.  OFF  ICE  SYMBOL 
ilf  applicable ) 

7*.  NAME  OF  MONITORING  ORGANIZATION 

Air  Force  Office  of  Scientific  Research 

6c.  AO  OR  ESS  <ca>.  State  and  /.IP  Code) 

Texas  A&M  University 

College  Station,  TX  77843 

Ba  NAME  CF  FUNDING-SPONSORING 
ORGANIZATION 

Bb.  OFFICE  SYMBOL 
(If  applicable) 

Be  ADDRESS  City,  State  and  /.If*  Code) 


7b.  ADDRESS  (City.  State  and  ZIP  Code) 

Bolling  AFB 

Washington,  D.C.  20332 


Grant  No.  AFOSR-84-0067 


10  SOURCE  OF  FUNDING  NOS. 


PROGRAM 
E  LEMENT  NO. 


11  TiTlE  ilnctude  Security  Classification) 

A  Damage  Model  For  Continuous  Fiber  Composites:  Part  I 


12.  PERSONAL  AUTHOR(S) 

D.H.  Allen,  S.E.  Groves,  R.A.  Schaper 


PROJECT 

TASK 

NO. 

NO. 

13a.  TYPE  OF  REPORT 

13b.  TIME  COVERED 

14.  DATE  OF  REPORT  iVr  ,  Mo.,  Day) 

Interim 

FROM  NA  TO 

FEB  1985 

15.  PAGE  COUNT 


17 

COS AT  1 

CODES 

FIELD 

GROUP 

SUB  GR 

18  SUBJECT  TERMS  i Continue  on  reverse  if  necessary  and  identify  t>y  block  number) 

composites  failure 

damage  internal  state  variables 

laminate  analysis  plasticity  _ 


19.  ABSTRACT  iContinue  on  reverse  if  necessary  and  identify  by  block  number ) 

A  continuum  mechanics  approach  is  utilized  herein  to  develop  a  model  for  predicting 
the  thermomechanical  constitution  of  continuous  fiber  composites  subjected  to  both  mono¬ 
tonic  and  cyclic  fatigue  loading.  In  this  model  the  damage  is  characterized  by  a  set 
of  second  order  tensor  valued  internal  state  variables  representing  locally  averaged 
measures  of  specific  damage  states  such  as  matrix  cracks,  fiber-matrix  debonding, 
interlaminar  cracking,  or  any  other  damage  state.  Globally  averaged  history  dependent 
constitutive  equations  are  posed  utilizing  constraints  imposed  from  thermodynamics  with 
internal  state  variables  as  well  as  fracture  mechanics. 

In  Part  I  the  thermodynamics  with  internal  state  variables  is  constructed  and  it  is 
shown  that  suitable  definitions  of  the  locally  averaged  field  variables  will  lead  to 
equivalent  thermodynamic  constraints  on  a  scale  assumed  to  be  large  compared  to  the  scale 
of  the  damage.  Based  on  this  result  the  Helmholtz  free  energy  is  then  expanded  in  a 


20  DISTRIBUTION  A  V  Al  LA0I  LI  T  Y  OF  ABSTRACT 

UNCLASSIF  IE  D  UNLIMITE  D  3  SAME  AS  RPT  _i  OTIC  USERS  D 


22a  NAME  OF  RESPONSIBLE  i  N  L  i  V  I D  U  A  L 

Ma  j  .  D.  Glasgow 


21  ABSTRACT  SECURITY  CLASSIFICATION 

Unci ass if ied 


22b  TELEPHONE  NUMBER 
‘Include  l  rra  Code  t 


(202)  767-4937 


T-  ■  -W  .  -V  1  -V  %  v  -.  V  '  *-v  '  “  '  .  “-Is*  •  -  -  -  v  \  ,  -  „  -  .  -  .  •.  '.  »  .  ^T*  fw  *  .  •  :  •  .  _  ■'. 


A  DAMAGE  MODEL  FOR  CONTINUOUS  FIBER  COMPOSITES 
PART  I:  Theoretical  Development 

by 

D.H.  Allen 
S.E.  Groves 

Aerospace  Engineering  Department 
and 

R.A.  Schapery 

Civil  and  Aerospace  Engineering  Departments 
Texas  A&M  University 
College  Station,  Texas  778^3 


MM-5023-84-1 7 


August  1985 
Revised  February  1985 


.V-V-V. 


A  DAMAGE  MODEL  FOR  CONTINUOUS  FIBER  COMPOSITES 

PART  I:  Theoretical  Development 

by 

D.H.  Allen 
S.E.  Groves 
R.A.  Schapery 

ABSTRACT 

A  continuum  mechanics  approach  is  utilized  herein  to  develop  a 
model  for  predicting  the  thermomechanieal  constitution  of  continuous 
fiber  composites  subjected  to  both  monotonic  and  cyclic  fatigue  loading. 
In  this  model  the  damage  is  characterized  by  a  set  of  second  order 
tensor  valued  internal  state  variables  representing  locally  averaged 
measures  of  specific  damage  states  such  as  matrix  cracks,  fiber-matrix 
debonding,  interlaminar  cracking,  or  any  other  damage  state.  Globally 
averaged  history  dependent  constitutive  equations  are  posed  utilizing 
constraints  imposed  from  thermodynamics  with  internal  state  variables  as 
well  as  fracture  mechanics. 

In  Part  I  the  thermodynamics  with  internal  state  variables  is 
constructed  and  it  is  shown  that  suitable  definitions  of  the  locally 
averaged  field  variables  will  lead  to  equivalent  thermodynamic 
constraints  on  a  scale  assumed  to  be  large  compared  to  the  scale  of  the 
damage.  Based  on  this  result  the  Helmholtz  free  energy  is  then  expanded 
in  a  Taylor  series  expansion  in  terms  of  strain,  temperature,  and  damage 
to  obtain  the  stress-strain  relation  for  composites  with  internal  state 
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variables  representing  damage.  Finally,  an  internal  state  variable 
growth  law  is  proposed  for  matrix  cracking. 

In  Part  II  the  resulting  three  dimensional  tensor  equations  are 
simplified  using  material  symmetry  constraints  and  are  written  in 
engineering  notation.  The  resulting  constitutive  model  is  then  cast 
into  laminate  equations  and  an  example  problem  is  solved. 

It  is  concluded  that  although  the  model  requires  further 
development  and  extensive  experimental  verif icatiion  it  may  be  a  useful 
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tool  in  characterizing  the  thermomechanical  constitutive  behavior  of 
continuous  fiber  composites  with  damage. 


INTRODUCTION 


Ultimate  failure  of  continuous  fiber  composite  structural 
components  is  preceded  by  a  sequence  of  microstructural  and 
macrostructural  events  such  as  microvoid  growth,  matrix  cracking,  fiber- 
matrix  debonding,  interlaminar  cracking,  edge  delamination,  and  fiber 
fracture  which  are  all  loosely  termed  damage.  Considerable  experimental 


research  has  been  performed  in  the  last  decade  detailing  the  growth  of 
damage  under  both  monotonic  and  cyclic  loading  conditions  [1-7].  The 
significance  of  this  damage  lies  in  the  fact  that  numerous  global 


material  properties  such  as  stiffness  and  residual  strength  may  be 
substantially  altered  during  the  life  of  the  component,  as  shown  in  Fig. 
1  [8].  It  has  been  found  that  the  first  phase  of  fatigue  is  typified  by 
development  of  a  characteristic  damage  state  (CDS)  [9]  which  is  composed 
primarily  of  matrix  cracking  in  off-axis  plies.  During  the  second 
phase  of  damage  development  the  CDS  contributes  to  fiber-matrix 
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debonding,  delamination,  and  fiber  microbuckling.  These  phenomena  in 
turn  contribute  to  a  tertiary  damage  phase  in  which  edge  delamination 
and  fiber  fracture  lead  to  ultimate  failure  of  the  specimen  [6]. 

Analytical  modeling  of  the  damage  state  appears  to  be  only 
recently  studied.  The  earliest  attempts  fall  under  the  general  heading 
of  laminate  analysis,  in  which  various  empirical  schemes  have  been 
developed  to  discount  ply  properties  in  the  presence  of  damage  [10-12]. 
Axial  stiffness  reduction  and  stress  distribution  in  the  CDS  have  also 
been  predicted  using  a  one-dimensional  shear  lag  concept  [5].  Several 
researchers  have  obtained  solutions  for  effective  moduli  of  elastic 
bodies  with  distributed  cracks  [13-17].  In  the  case  where  cracks  are 
either  randomly  distributed  or  oriented  the  effect  of  total  crack 
surface  area  is  found  to  cause  a  first  order  effect  on  the  stiffness 
[13-16]. 

Fracture  based  concepts  have  recently  been  utilized  to  model 
damage  development  [18-21].  Although  the  first  of  these  studies  [18] 
contains  a  general  theory  which  may  be  applied  to  fibrous  composites,  it 
has  so  far  only  been  utilized  for  quasi-isotropic  random  particulate 
composites  such  as  solid  rocket  propellant,  [9]  and  as  such  has  not  been 
applied  to  continuous  fiber  composites.  The  theory  in  the  latter  two 
[20,21]  has  not  been  utilized  to  predict  reduction  of  off-axis  stiffness 
components.  Kachanov's  modulus  reduction  technique  [22]  has  also  been 
utilized  for  fibrous  composites  [23]  and  although  promising  results  were 
obtained,  the  model  was  constructed  in  uniaxial  form  only. 

A  complex  interactive  experiment  and  analysis  model  (called  a 
mechanistic  model)  has  been  proposed  [8].  The  approach  used  therein  is 
fundamentally  quite  different  from  that  developed  in  this  paper. 


Furthermore,  the  mechanistic  model  requires  numerous  experimental 


results  for  each  geometric  layup  in  order  to  determine  which  damage  mode 
predominates. 

Finally,  extended  forms  of  Miner's  rule  [24]  have  been  proposed 
[18,25].  However,  they  are  based  on  simplified  microphysical  models  at 
this  time. 

The  concept  of  damage  as  an  internal  state  variable  [26]  has  been 
previously  utilized  in  continuum  mechanics/thermodynamics  based  theories 
for  crystalline  and/or  brittle  materials  [27-34],  as  well  as  for 
nonlinear  viscoelastic  materials  [18].  A  study  has  been  made  of  the 

t 

effect  of  vector-valued  damage  parameters  on  various  compliance  terms 
[35],  and  this  methodology  is  currently  undergoing  further  development 
[36,  37]. 

The  research  reviewed  above  indicates  that  although  important 
progress  has  been  made  in  characterizing  damage  in  fibrous  composites, 
substantial  and  continued  research  is  warranted  before  several  issues 
can  be  resolved. 

In  this  paper  an  attempt  will  be  made  to  assemble  many  of  the 
concepts  embodied  in  the  research  efforts  mentioned  briefly  above  and  to 
utilize  these  concepts  to  develop  a  model  for  damage  in  continuous  fiber 
composites  which  is  rigorously  based  in  continuum 
mechanics/thermodynamics  and  is  generic  with  regard  to  material  type, 
load  spectrum,  and  specimen  geometry. 
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CHARACTERIZATION  OF  DAMAGE  AS  A 


SET  OF  INTERNAL  STATE  VARIABLES 

Consider  an  initially  unloaded  and  undamaged  continuous  fiber 
composite  structural  component  as  shown  in  Fig.  2a,  where  undamaged  is 
defined  here  to  mean  that  the  body  may  be  considered  to  be  continuous 
(without  voids)  on  a  scale  several  orders  of  magnitude  smaller  than  the 
smallest  external  dimension  of  the  component.  Although  voids  may  exist 
in  the  initial  state,  their  total  surface  area  is  assumed  to  be  small 
compared  to  the  external  surface  area  of  the  component.  Under  this 
assumption  the  body  is  assumed  to  be  simply  connected  and  we  call  the 
initial  bounding  surface  the  external  boundary.  Although  the  component 
is  undamaged,  there  may  exist  local  heterogeneity  due  to  such  causes  as 
processing  inhomogeneities  and  second  phase  materials  including  fibers 
and  matrix  tougheners.  In  addition,  the  body  may  be  subjected  to  some 
residual  stress  state  due  to  prior  loading,  cool  down,  etc. 

Now  suppose  that  the  component  is  subjected  to  some  traction  or 
deformation  history,  as  shown  in  Fig.  2b.  The  specimen  will  undergo  a 
thermodynamic  process  which  will  in  general  be  in  some  measure 
irreversible.  This  irreversibility  is  introduced  by  the  occurrence  of 
such  phenomena  as  material  inelasticity  (even  in  the  absence  of  damage) 
fracture  (both  micro-  and  macroscale),  friction  (due  to  rubbing  of 
fractured  surfaces),  and  chemical  change.  While  all  of  these  phenomena 
can  and  do  commonly  occur  in  continuous  fiber  composites,  in  the  present 
research  it  will  be  assumed  that  fracture  is  the  only  irreversible 
phenomenon  of  significance.  Thus,  all  fracture  events  will  be  termed 
damage.  Due  to  these  fracture  events,  the  body  will  necessarily  become 
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multiply  connected,  and  all  newly  created  surfaces  not  intersecting  the 

external  boundary  will  be  termed  internal  boundaries.  Because  of  the 

above  assumptions  the  model  will  be  limited  to  polymeric  matrix 

composites  at  temperatures  well  below  the  glass  transition  temperature 

T  ,  where  viscoelasticity  in  matrix  materials  is  small.  Metal  matrix 
S 

composites  will  be  excluded  due  to  complex  post-yield  behavior  of  the 
matrix. 

While  fracture  classically  involves  changes  in  the  boundary 
conditions  governing  a  complex  field  problem,  it  is  hypothesized  that 
one  may  neglect  boundary  condition  changes  caused  by*  creation  and 
alteration  of  both  internal  and  external  surfaces  created  during 
fracture  as  long  as  the  resulting  damage  in  the  specimen  is 

statistically  homogeneous  on  a  scale  which  is  small  compared  to  the 

scale  of  the  body  of  interest.  This  implies  that  the  total  newly 

created  surface  area  (which  includes  internal  surfaces)  may  be  of  the 
same  order  of  magnitude  as  the  original  external  surface  area.  Under 
the  condition  of  small  scale  statistical  homogeneity  all  continuum  based 
conservation  laws  are  assumed  to  be  valid  on  a  global  scale  in  the  sense 
that  all  changes  in  the  continuum  problem  resulting  from  internal  damage 
are  reflected  only  through  alterations  in  constitutive  behavior. 
Typical  microstructural  events  which  qualify  as  damage  are  therefore 
matrix  cracking  in  lamina,  fiber/matrix  debonding,  and  localized 
interlaminar  delamination.  Large  scale  changes  in  the  external  surface 
such  as  edge  delaminations,  however,  are  treated  as  boundary  effects 

which  must  be  reflected  in  conservation  laws  via  changes  in  the  external 
boundary  conditions  rather  than  in  constitutive  equations  [39]. 
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Although  the  damage  process  actually  involves  the  conversion  of 
strain  energy  to  surface  energy,  the  fact  that  the  damage  is  reflected 
in  the  local  constitutive  equations  requires  that  it  be  treated  as  a  set 
of  energy  dissipative  local  state  variables  which  are  not  discernible  on 
the  external  boundary.  Therefore,  since  the  damage  can  be  determined 

only  through  a  precise  knowledge  of  the  entire  history  of  observable 

inputs,  it  is  characterized  as  a  set  of  internal  state  variables.  This 
concept  will  be  further  developed  in  the  next  section. 

THERMODYNAMICS  OF  MEDIA 
WITH  DAMAGE 

We  now  proceed  to  construct  a  concise  model  of  the  continuum  with 

damage.  To  do  this,  consider  once  again  the  structural  component, 

denoted  B  in  Fig.  3a.  The  body  B  is  assumed  to  be  of  the  scale  of  some 

appropriate  boundary  value  problem  of  interest.  Now  consider  some  local 

element  labelled  L  and  with  external  surfaces  arbitrarily  chosen 

normal  to  a  set  of  Cartesian  coordinate  axes  (x^,  x^,  x^),  as  shown  in 

Fig.  3b.  The  element  L  is  extracted  from  B  and  the  newly  created 

surfaces  are  subjected  to  appropriate  boundary  conditions  so  that  the 

element  response  is  identical  to  that  when  it  is  in  B.  Internal 

surfaces  caused  by  fracture  are  labelled  such  that  the  intersection 

of  S1  and  S,,  is  a  null  set  and  S  =  +  S^.  Furthermore,  the  volume  of 

the  element  is  defined  to  be  V.  ,  which  includes  the  volume  of  voids  V  . 

L  v 

The  scale  of  L  is  chosen  so  that  its  dimensions  are  small  compared  to 
the  dimensions  of  B,  but  at  the  same  time,  the  dimensions  of  L  are  large 
enough  to  guarantee  statistical  homogeneity  of  the  material 


inhomogeneities  and  defects  in  L  even  though  the  total  surface  area  of 
defects  may  be  of  the  same  order  of  magnitude  as  S1  [38].  Suppose 
further  that  in  the  absence  of  defects  or  at  constant  damage  state 
the  material  behavior  is  linear  thermoelastic,  thus  specifically 
excluding  the  effects  of  crack  face  rubbing. 

The  following  notation  is  adopted.  Quantities  without  capitalized 
subscripts  denote  pointwise  quantities.  Those  with  subscripts  L  denote 
quantities  which  are  averaged  over  the  local  element  L.  Finally,  the 
subscript  E  denotes  linear  thermoelastic  properties. 

Review  of  Thermodynamic  Constraints  on  Linear  Thermoelastic  Media 

Under  the  above  conditions  the  pointwise  Helmholtz  free  energy  per 
unit  mass  h  of  the  undamageu  linear  elastic  medium  may  be  expressed  as  a 
second  order  expansion  in  terms  of  strain  and  temperature  change  AT 
as  follows  [40]: 

h  =  u  -  Ts  =  h£  (eijtT)= 

A  +  3.  e.  .  +  1/2  C .  . .  .  e .  . e.  .  +  DAT  +  E. .e. .AT  +  1/2  FAT2  (1 ) 

ij  lj  ljkl  lj  kl  lj  lj 

where  u  and  s  are  the  internal  energy  and  entropy  per  unit  mass, 
respectively,  and  A,  B.  ,  C.^^,  D»  E  and  F  are  material  parameters 
which  are  independent  of  strain  and  temperature  and  AT  =  T-T_,  where  TD 

A  A 

is  the  reference  temperature  at  which  no  deformation  is  observed  at  zero 

load.  In  addition,  we  assume  here  that  all  motions  produce  small 

deformations  so  that  e..  is  the  infinitesimal  strain  tensor. 

ij 

Furthermore,  inertial  effects  are  assumed  to  be  negligible. 

Pointwise  conservation  laws  appropriate  to  the  body  are  as 
follows: 


1)  conservation  of  linear  momentum 


(2) 


o  .  . ,  .=0 
Ji’  J 


where  is  the  work  conjugate  stress  tensor  to  the  strain  tensor  e . ^ 
and  body  forces  are  assumed  to  be  negligible; 

2)  conservation  of  angular  momentum  (assuming  body  moments  may  be 
neglected) 

°u  '  °ji  •  <3> 

3)  balance  of  energy 


pu  -  o .  .e .  .  +  q .  .  =  pr 
ij  ij  J,J 


(4) 


where  p  is  the  mass  density,  q^  are  the  components  of  the  heat  flux 

vector,  and  r  is  the  heat  source  term;  in  addition,  dots  denote  time 

differentiation  and,  .  ~3/3x.; 

.J  y 

4)  the  second  law  of  thermodynamics 

ps  -  £1  +  Ak  .  1  0  ,  1  (5) 

T  '-x  J  *  3 
Furthermore 

e,  ,  :  1  (u,  .  +  u  .  .  )  , 

IJ  2  J’1  (6) 

The  above  set  of  equations  may  be  cast  with  appropriate  boundary 
conditions  so  that  constraints  imposed  by  the  second  law  of 

thermodynamics  will  result  in  [40] 

-3h„ 


s  =  sE  = 


3T 


E  =  -D-E  .e. .-FAT 
ij  iJ 


(7) 


a. .  =  a 
IJ  E. 


3hF 

=  p  E 


IJ 


3e 


ij 


0  (Bij  +  Cij WElct  +  Eij  4T> 


(8) 


where  pB . ^  are  components  of  residual  stresses  at  "he  reference 


temperature  at  which  T  =  0, 
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(9) 


'  "CT  '-K  \  *  K  i  ^  " 


i  k  wwfl  vip. 


=  -k 


ij°j 


+  0  ( 


12s 


where 


gj 


(10) 


Thermodynamic  Constraints  with  Local  Damage 

Now  consider  the  local  element  shown  in  Fig.  3b  with  traction 
boundary  conditions  on  the  external  surface  S1  .  In*  addition,  the 
interior  of  L  is  assumed  to  be  composed  entirely  of  linear  elastic 
material  and  voids.  It  is  our  aim  to  construct  locally  averaged  field 
equations  which  are  similar  in  form  to  the  pointwise  field  equations 
discussed  above.  In  performing  this  averaging  process  the  pointwise 
Helmholtz  free  energy  described  in  equation  (1)  will  undergo  a  natural 
modification  to  include  the  energy  conversion  due  to  crack  formation. 

In  order  to  construct  locally  averaged  equations  first  define  the 
locally  averaged  stress  tensor: 


(11) 


Integrating  the  balance  of  energy  (4)  over  the  local  volume  and  dividing 
through  by  the  local  volume  results  in 
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In  order  to  simplify  the  above  define  the  locally  by  averaged  density  p 

Li 

such  that 


Furthermore,  define  the  local  internal  energy  rate  per  unit  mass  it  such 

Li 

that 


which  can  be  constructed  as  long  as  pL  t  0. 

Now  consider  the  second  term  in  equation  (12).  Recall  that  since 
o . j  is  a  symmetric  tensor 


a.  .c.  . 
ij  ij 


2 


a.,  (u.  .  +  u .  . )  =  a . .u .  . 

ij  i.J  J.i  ij  i>J 


Thus,  using  the  divergence  theorem  gives 


..  .  t/'M. 


cLs~inrf 


Substituting  equations  (14),  (20),  (21),  and  (22)  into  , equation  (12) 

yields  the  following  locally  averaged  balance  of  energy: 


PLUL+  plA  "  °L.  A.  .  +  qL 


=  PtrT 


1J  1J  J.J 


which  can  be  seen  to  be  equivalent  in  form  to  energy  balance  law  (4), 


irAWi 


Define  also 


We  now  define  the  effective  internal  energy  uL  such  that 


UL  -=  UL  +  UL 


Substitution  of  (24)  into  (23)  results  in 


PLUL  ~  %  ,eL..  +  %.  .  =  PLrL 
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In  order  to  construct  a  similar  statement  for  entropy  production 
ality  (5),  first  multiply  through  by  T  and  then  integrate  over  the 


inequa 


local  volume  V  and  divide  by  this  quantity  to  obtain 

Li 


+  i  /T(^. 


,  dV  >  0 
►  J  ~~ 


Now  define 


'V 


so  that  substitution  of  definitions  (22),  (27)  and  (28)  into  (26)  will 
result  in 


PLSLTL  -  pLrL  + 


dV  >  0 


Now  note  that  the  last  term  in  (29)  may  be  written  as  follows  using  the 
product  rule: 


[t(^)  ,  dv  =  ^  Tq,  =dv  -  —  f  qj 


Thus,  define 


T 


L, 


j 


(31) 


Therefore,  substituting  definitions  (21)  and  (31)  into  (30)  and  this 

result  into  (29)  gives,  after  dividing  through  by  T  : 

L 


similar  to  pointwise  equations  (2),  and  the  conservation  of  angular 
momentum  may  also  be  obtained 


similar  to  equations  (3).  Thus,  it  is  assumed  that  no  body  moments  are 


introduced  via  material  inhomogeneity  or  other  sources. 

Locally  averaged  field  equations  (25),  (32), (33)  and  (3*0  have  now 
been  constructed  which  are  similar  to  pointwise  field  equations  (2) 
through  (5).  On  the  basis  of  this  similarity  we  now  define  the  locally 
averaged  Helmoltz  free  energy  [18,  38]: 


tlsl 


UL 


tlsl+  ul 


=  h„ 


where  it  can  be  seen  from  definitions  (14)  and  (28)  that  h„  is  the 
locally  averaged  elastic  Helmholtz  free  energy  for  which  residual 
damage  is  assumed  to  be  small  ,  given  by  equation  (1)  to  be 


if. 


(A  +  B  .  c.  .  +  1  C.  ..  »e.  .  e. 
v  ij  ij  T  13  k  l 


HAT  +  E. ,e. . AT+1  FAT  ) 
ij  ij  ~ 


Utilizing  definitions  (18)  and  (27)  we  may  further  define 


1 


no  sum  on  i , j , 


C,  eT  CT  -  1  /  C...,  e,-4  CV1  dv  »  no  sum  on  i>  i  •  k>  x> 

L.  .  L  v  I  KI 


dV 


•V  •  YV -V LWV.M 


DI  =  _L  f 

V  T  J 
L  L-'V 

L 

EL. .  CL. .  =  -- 
13  1J  V. 


DATdV 


v/ 

LL-'V, 


E..  e..  TdV  ,  and  sum  on  i,j,  and 

13  13 


F  =  1  f  FT2dV 

W  Jvl 

Note  that  when  the  actual  strain  state  e  is  spatially  invariable  and 

S_  =0,  the  above  equations  result  in  B.  =B. . ,  C.  =  C. ,  and  E. 

2  Lij  iJ  Lijkl  1Jkl  L] 

=  E_ .  Substituting  definitions  (37)  through  (42)  into  local  elastic 


free  energy  equation  (36)  results  in 


X  '  ^  +  + 1  +  +  ♦ 7  v,2  •  <“> 


The  similarity  between  the  pointwise  and  local  field  equations 


leads  to  the  conclusion  that 


L  9T 


'•  a 

'v'yv 

vv.V- 


°L. .  PL  3e. 


pL  <bl  +  Cl  eL  +  EL  .V  +  PL  3eT 

L  Lij  ijkX.  Lk l  ij  I 


,  =  -kL  .8L.  +  °  (l8L.!  }  ’ 

i  13  3  J 


where 


dV 


no  sum  on  i»j» 


l*»  - t»  i‘a  '.V  V.  r'm.  tV 
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r*r\  /  "  t  r;  v* 


> -\.*h 
tAjft 

tv* 

V-V. 

*■.*■*.•' 

w'.*-'.*i 


and 


Li  =  TL,j 


(48) 


*V  v 

;v> 


r.-TJ? 


Note  the  similarity  between  equations  (7)  through  (10)  and  (44)  through 
(48),  respectively. 

Equations  (45)  will  serve  as  the  basis  for  thermomechanical 
stress-strain  relations  in  damaged  fibrous  composites.  All  damage  will 
be  reflected  through  the  local  energy  due  to  cracking  u^.  This  term 
will  be  modelled  with  internal  state  variables  characterizing  the 
various  damage  modes. 


I  ■  •  V 


Description  of  the  Internal  State 

In  order  to  describe  the  internal  state,  we  first  consider  the 
kinematics  of  a  typical  point  0  with  neighboring  points  A  and  B,  as 
shown  in  Fig.  4.  Before  deformation  lines  0A  and  OB  are  orthogonal,  as 
shown  in  (a).  After  deformation  we  imagine  that  lines  joining  O’, A’, 
and  B'  are  as  shown  in  (b),  and  just  at  the  instant  that  deformation  is 
completed,  a  crack  forms  normal  to  the  plane  of  AOB  through  point  O',  as 
shown  in  (c).  Furthermore,  point  O'  becomes  two  material  points  O'  and 
0”  on  opposite  crack  faces  and  points  A'  and  B'  deform  further  to  points 
A"  and  B" .  It  is  assumed  that  all  displacements,  including  crack 
opening,  are  infinitesimal,  so  that  an  observer  at  an  appropriate 
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observation  distance  from  point  0  "sees"  only  the  deformation  A"  O'  B". 


The  strain  associated  with  this  deformation  is  appropriately  called  an 
observable  state  variable.  However,  the  strain  of  interest  is 
associated  with  A"0"B".  Therefore,  it  is  essential  to  construct  an 

internal  state  variable  which  will  relate  these  two  strain  descriptions. 

Y  Q  _*.Q 

We  therefore  construct  the  vectors  u  connecting  O'  and  0"  and  n 

describing  the  normal  to  the  crack  face  at  O',  as  shown  in  Fig.  5.  It 

jn 

should  be  noted  that  u  can  be  used  to  construct  a  pseudo-strain 

representing  the  difference  in  rotation  and  extension  of  lines  A"0'B" 

and  A"0"B" . 

Now  recall  that  the  energy  released  during  cracking  is  given  by 
equation  (17).  Since  the  body  is  elastic,  we  assume  that  this  process 
is  reversible  and  that  tractions  T^  can  be  applied  at  point  O'  which 
will  close  the  crack: 


1 


•c 

UT  =  „ 

L  pT  V 


L  L 


(49) 


Using  Cauchy's  formula  the  above  can  be  placed  in  a  form  similar  to 
(17): 


1 


P  V 
y  L  L 


(50) 


where  the  superscripts  denote  quantities  associated  with  the  actual 
crack  geometry.  Although  these  quantities  do  not  necessarily  coincide 
with  the  terms  in  the  integrand  of  (17),  their  surface  integration  will 


( 


'  *  * 


< 
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result  in  precisely  the  same  energy  release  rate  due  to  the 

reversibility  of  the  process. 

-h.C 

Guided  by  the  fact  that  u  and  n  describe  the  kinematics  of  the 
cracking  process  at  point  0,  we  now  define  the  following  second  order 
tensor  valued  internal  state  variable: 


“u  3  \  »' 


Note  that  when  the  crack  normal  is  time  independent  the  above  may  be 
time  differentiated  to  give 


a  .  .  a  .  .  dS 
iJ  ij 


Note  that  the  components  of  u  and  n  can  be  recovered  from  (51) 
by  using  single  row  and  column  multiplication  of  : 


2  c  c 
u.  =  u.  r.  u_.  n 


(no  sum  on  i) , 


2  =  uc  n  u.  n. 
"i  i  J  1  J 


(no  sum  on  j ) . 


Therefore,  the  normal  and  shear  modes  of  crack  displacement  can  be 

recovered  from  a. .. 

ij 

Now  suppose  that  uc  is  subdivided  into  integrals  over  a  finite 


number  of  internal  surface  areas  S  '  of  fundamentally  different  nature: 


such  that 


and  the  integer  N  represents  the  number  of  damage  modes,  to  be  described 
below. 

Now  define  the  locally  averaged  internal  state  variable  for 

the  n  th  damage  mode  as  follows: 


Furthermore,  if  we  define  the  average  crack  closure  stress 
for  the  n  th  damage  mode  such  that 


no  sum  on  i ,  j ,  n , 


it  follows  that,  from  (55),  (57)  and  (58) 


where  we  have  assumed  that  repeated  indices  n  imply  summation  over  the 


range  N.  It  is  clear  from  the  above  discussion  that  the  value  of  N  must 
be  sufficiently  large  to  recover  the  essential  physics  of  the  damage 


process.  In  a  mathematical  sense,  this  implies  that,  whereas  the 


mapping  from  a. .  to  a?  is  unique,  the  inverse  should  also  be  true  in 
1 J  L* .  . 

1J 


an  approximate  sense.  Note  also  that  both  u^  and  n^  in .equations  (57) 
will  be  affected  by  crack  interaction  in  the  local  volume. 


Now  consider  equation  (59)  in  further  detail.  The  kinetic 
cn 


quantities  o 


may  be  interpreted  as  generalized  forces  which  are 


ij 


energy  conjugates  to  the  kinematic  strain-like  internal  state  variables 


a, 


n 


V 


We  infer  from  this  that  there  exists  a  constitutive  relation,  between 
these  variables  of  the  form 


cn 


=  oT 


cn 


<eT  .  T  ,  a,"  ) 


ij 


“kl 


kl 


Therefore,  substituting  (60)  into  (59)  and  integrating  in  time 
will  give 


uC  (t,  )  =  f  u7C  (t>dt  "  UT  (eL.,  (tl>’  tL  (i]). 


(t,))  , 


It  is  now  proposed  that  u  be  expanded  in  a  Taylor  series  which  is 
second  order  in  each  of  the  arguments  in  equation  (61)  as  follows: 
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where  all  terms  are  at  least  linear  in  due  to  the  fact  that 

ij 

depends  explicitly  on  damage,  and  aT^  =  TL~TR.  Thus,  substituting  (43) 
and  (62)  into  equations  (45)  and  neglecting  higher  order  terms  yields: 


<*L.  .  =  B  +  Et  ATt  +  CT  e  +  I.  ..  . 

!!  Lij  Lij  L  Lijkl  Lk]^  ijkl  Lkl 
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Equations  (63)  may  be  written  in  the  following  convenient  form 


where 


+ 


n 

I 

ijkl 


n 


a 


kl 


is  the  residual  stress  in  the  absence  of  strain  and  temperature  change, 
which  may  be  induced  by  damage; 


c'  h  C  +  Mn  a? 

Lijkl  Lijkl  ijklmn  Lmn 


is  called  the  effective  modulus  tensor,  which  is  degraded  due  to  damage; 
and 


is  called  the  thermal  strain,  which  also  may  be  affected  by  damage. 
Note  that  the  only  second  order  terms  are  due  to  damage  induced 
stiffness  reduction  and  thermal  strains. 

Equations  (6^)  through  (67)  are  the  completed  description  of  the 
stress-strain  relationship.  Note  that  these  equations  reduce  to  the 
standard  linear  thermoelastic  equations  in  the  absence  of  damage  (a 
=0).  Furthermore,  it  is  notable  that  these  equations  reduce  to 


Kachanov's  model  [22]  in  uniaxial  form. 


Damage  Growth  Laws 


The  model  is  completed  with  the  construction  of  the  damage  growth 
laws,  which  may  be  described  in  the  following  differential  equation 
form: 


;n  =  ^  (C  T  ' 

ij  ULkl’  TL  hk] 


or  equivalently,  when  1  '  are  analytic  in  time, 

ij 

“E  ,<V  ■  f  '  Xi  ■  . 

ij  J 

—  CO 

where  the  dot  denotes  time  differentiation.  Although  the  above 
equations  are  called  "growth"  laws  they  have  the  more  general  capability 
to  model  such  phenomena  as  healing. 

The  precise  nature  of  equations  (69)  is  determinable  only  through 
a  concise  experimental  program  coupled  with  an  understanding  of  the 
micromechanics  of  the  medium.  Indeed,  these  growth  laws  constitute  the 
single  most  complex  link  in  the  model  development. 

In  this  section  an  example  of  a  first  generation  growth  law  will 
be  constructed  for  predicting  damage  up  to  the  CDS  in  continuous  fiber 
composites.  Experimental  evidence  suggests  that  matrix  cracks 

predominate  prior  to  development  of  the  CDS  [  *4  — 6] .  Guided  by  this 

1 

observation,  a  single  damage  tensor  is  considered  herein:  a. 

Lj  .  . 
il 

representing  matrix  cracking. 

In  order  to  completely  define  equations  (69),  it  is  necessary  to 
construct  indicators  of  both  the  magnitude  and  direction  of  the  damage 


tensor.  In  this  first  generation  model  it  is  assumed  that  the  direction 


of  the  damage  tensor  is  known  a  priori  and  does  not  vary  as  the  damage 


state  changes.  Specifically,  in  a  typical  laminate,  it  is  assumed  that, 


in  accordance  with  equation  (57),  the  locally  averaged  resultants  of  u 


and  n  are  normal  to  the  fiber  direction  in  each  ply,  as  shown  in  Fig. 


6.  Thus,  for  example,  in  a  0°  ply  oJ  0,  and  all  other  components 


are  zero,  whereas  in  a  90°  ply,  a}  0  and  all  other  components  are 

L11 

zero  in  global  coordinates. 


Under  the  above  assumptions,  the  magnitude  of  the  damage  tensor  is 


the  sole  repository  for  history  dependence  in  each  ply.  Experimental 


evidence  indicates  that  for  matrix  cracking  in  randomly  oriented 


particulate  composites  [*41]  and  matrix  cracks  in  fibrous  composites 


[20,21]  the  growth  of  damage  surface  area  is  related  to  the  energy 


release  rate  G  by 


da  n 

_  a  G 

dN 


where  N  represents  the  number  of  cycles  in  a  fatigue  test,  and  n  is  some 


material  parameter.  Guided  by  these  results,  a  similar  law  is 


constructed  here.  To  accomplish  this,  first  multiply  both  sides  of  (70)  by 


dN/dt  to  obtain  for  the  90°  ply,  for  example, 


cut,  dN 
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Assuming  that  the  energy  release  rate  is  essentially  mode  I  and 


therefore  depends  on  the  maximum  normal  strain,  the  damage  growth  law 
for  transverse  cracking  is  thus  hypothesized  to  be  of  the  form 


K 


c  <  e  <  e  ,  and 
n  .  n  n 
min  max 


o  if  e, 


<  e  or 

n  , 
min 


e  >  e 

n  n 

max , 


(72) 


where  is  the  local  normal  strain  component  which  is  normal  to  the 

fibers.  Furthermore,  e  is  the  value  of  e  at  which  transverse 

nmin  n 

cracking  initiates,  and  e  is  the  corresponding  value  at  which 

max 

transverse  cracking  saturates.  K  and  n  are  experimentally  determined 
material  parameters  which  may  depend  on  the  initial  damage  state  or  on 
history  dependent  damage  such  as  matrix  microvoids.  The  use  of 
presupposes  that  the  fracture  mode  is  predominantly  mode  I  in  nature, 
which  may  not  be  the  case  in  some  complex  layups.  In  these  cases,  mode 
II  and  mode  III  terms  may  be  required.  Note  that  all  components  of  a 

L  .  . 

ij 

are  zero  except  a.  ,  which  is  nonzero  in  the  local  coordinate  system 
2  2 

wherein  the  fibers  are  aligned  parallel  to  the  x^  axis. 

The  time  derivative  of  reflects  the  time  independent  nature  of 

the  damage  mode,  and  actually  represents  a  departure  from  the  form  of 

equations  (69)  in  that  explicit  strain  rate  dependence  is  now  reflected 

in  the  growth  law.  However,  damage  growth  law  (73)  need  not  satisfy  the 

principal  of  equipresence  [42]  since  it  is  not  an  equation  of  state 
[43]  . 


Futhermore,  a  study  of  thermodynamic  constraints  with  linear  strain  rate  independence 
in  equations  (69)  will  indicate  that  equations  (44)  through  (46)  remain  valid  [44,45]. 

Fig.  7  shows  a  typical  growth  history  for  a  specimen  subjected  to  monotonically 
increasing  deformation  u(L). 

Equation  (72)  completes  the  description  of  the  damage  model.  Integration  of 
these  equations  in  time  will  lead  to  current  values  of  the  damage  tensor  which  is 
input  to  constitutive  equations  (64)  through  (67).  It  should  be  pointed  out,  however, 
that  these  equations  may  be  extremely  nonlinear  and  as  such  must  in  many  cases  be 
integrated  numerically  with  stiff  integration  schemes  [48] . 


CONCLUSION 


Stres-strain  relations  have  been  developed  herein  which  account 
for  various  forms  of  damage  in  continuous  fiber  composites. 
Furthermore,  a  damage  growth  law  has  been  proposed  for  matrix  cracking 
in  fibrous  composites.  The  model  developed  herein  is  thus  a  complete 
description  necessary  to  characterize  the  thermomechanical  constitution 
of  a  fibrous  composite  (excluding  failure).  However,  the  actual  use  of 
this  model  is  complicated  by  the  requirement  for  numerous  experimentally 
determined  quanitites,  as  well  as  the  necessity  to  determine  locally 
based  observable  state  variables  by  analytic  methods.  The  construction 
of  these  parameters  constitutes  an  entire  separate  research  effort  which 
is  considered  in  Part  II. 
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constitutive  equations  in  the  standard  laminate  coordinate  system,  a  specified  constitutive 
model  is  developed  for  the  case  of  matrix  cracks  only.  The  potential  of  the  model  to 
predict  degraded  or  effective  stiffness  moduli  is  demonstrated  by  solving  the  problem  of 
transverse  matrix  cracks  in  the  9 (f  layer  of  [0/90]  and  [0/90„]s  laminates. 

To  solve  the  example  problems,  the  undamaged  modulii  are  determined  from  experimental 
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a  variety  of  matrix  crack  spacings  in  the  9CP  layers.  The  internal  state  variable  for  •/. 
transverse  matrix  cracking  is  related  to  the  strain  energy  release  rate  due  to  cracking  by 
utilizing  linear  elastic  fracture  mechanics.  The  values  of  effective  (damage  degraded) 
stiffness  predicted  by  the  constitutive  model  are  in  close  agreement  with  both  experiments 
and  finite  element  results.  The  close  agreement  obtained  in  these  example  problems,  whiliL 
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model  to  predict  degraded  stiffness. 
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A  DAMAGE  MODEL  FOR  CONTINUOUS  FIBER  COMPOSITES 
PART  II:  MODEL  APPLICATIONS 

by 

S.E.  Groves 
D.H.  Allen 
C.E.  Harris 
and 

R.A.  Schapery 

ABSTRACT 

A  continuum  mechanics  approach  is  utilized  herein  to  develop  a 
model  for  predicting  the  thermomechanical  constitution  of  continuous 
fiber  composites  subjected  to  both  monotonic  and  cyclic  fatigue  loading. 
In  this  model  the  damage  is  characterized  by  a  set  of  second  order 
tensor  valued  internal  state  variables  representing  locally  averaged 
measures  of  specific  damage  states  such  as  matrix  cracks,  fiber-matrix 
debonding,  interlaminar  cracking,  or  any  other  damage  state.  Locally 
averaged  history  dependent  constitutive  equations  are  posed  utilizing 
constraints  imposed  by  thermodynamics  with  internal  state  variables. 

In  Part  I  the  thermodynamics  with  internal  state  variables  was 
constructed.  It  was  shown  that  suitable  definitions  of  the  locally 
averaged  field  variables  led  to  equivalent  thermodynamic  constraints  on 
a  scale  assumed  to  be  large  compared  to  the  scale  of  the  damage.  Based 
on  this  result  the  Helmholtz  free  energy  was  expanded  in  a  Taylor  series 
in  terms  of  strain,  temperature,  and  damage  to  obtain  the  stress-strain 
relations  for  composites  with  internal  state  variables  representing 
damage.  Finally,  an  internal  state  variable  growth  law  was  proposed  for 
matrix  cracking. 
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In  this  paper,  the  three  dimensional  tensor  equations  from  Part  I 
[1]  are  simplified  using  symmetry  constraints.  After  introducing 
engineering  notation  and  expressing  the  constitutive  equations  in  the 
standard  laminate  coordinate  system,  a  specialized  constitutive  model  is 
developed  for  the  case  of  matrix  cracks  only.  The  potential  of  the 
model  to  predict  degraded  or  effective  stiffness  moduli  is  demonstrated 
by  solving  the  problem  of  transverse  matrix  cracks  in  the  90°  layer  of 
[0/90]  and  [0/90.,]  laminates. 

To  solve  the  example  problems,  the  undamaged  moduli  are 
determined  from  experimental  data.  The  damage  tensor  is  determined 
analytically  from  a  finite  element  analysis  assuming  a  variety  of  matrix 
crack  spacings  in  the  90°  layers.  The  internal  state  variable  for 
transverse  matrix  cracking  is  related  to  the  strain  energy  release  rate 
due  to  cracking  by  utilizing  linear  elastic  fracture  mechanics.  The 
values  of  effective  (damage  degraded)  stiffnesses  predicted  by  the 
constitutive  model  are  in  close  agreement  with  both  experimental  and 
finite  element  results.  The  close  agreement  obtained  in  these  example 
problems,  while  limited  to  transverse  matrix  cracks  only,  demonstrates 
the  potential  of  the  constitutive  model  to  predict  degraded  stiffnesses. 
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INTRODUCTION 


In  Part  I  [1]  it  was  hypothesized  that  damage  can  be  modeled  in 
continuous  fiber  composites  by  a  set  of  internal  state  variables  which 
represent  locally  averaged  measures  of  matrix  cracking,  interlaminar 
delar:  ination,  and  other  damage  mechanisms  on  a  scale  which  is  assumed  to 
be  small  compared  to  the  scale  of  the  boundary  value  problem  of 


2 


interest.  Continuum  mechanics  with  internal  state  variables  [1]  was 
then  utilized  to  construct  stress-strain  relations  in  which  all 
components  of  the  degraded  modulus  tensor  can  be  determined  for  a  given 
damage  state. 

The  purpose  of  this  paper  (Part  II)  is  to  demonstrate  how  the 
model  may  be  utilized  to  predict  the  stiffness  of  laminates  which  are 
subjected  to  known  damage  states.  This  procedure  is  illustrated  via 
specific  examples  in  which  there  is  a  single  damage  mode  consisting  of 
transverse  matrix  cracking.  It  is  shown  that  single  lamina  properties 
in  the  presence  of  damage  can  be  utilized  as  given  properties  to  obtain 

i 

favorable  comparisons  to  both  experimental  and  finite  element  results 
for  specific  laminates. 

The  model  application  is  accomplished  by  first  imposing  symmetry 
constraints,  performing  the  laminate  integration,  and  finally  reducing 
to  generalized  plane  stress. 

SIMPLIFICATION  OF  THE  MODEL 

We  now  consider  the  stress-strain  relation  described  in  equations 
(64)  through  (67)  of  Part  I  (see  Appendix  A).  For  the  examples  to  be 
considered  herein,  it  is  assumed  that  all  residual  stress  components  are 

P 

zero  (o  =0),  and  that  there  are  no  thermal  transients  (AT  =0). 

L.  .  L 

ij 

Reduction  to  Single-Index  Notation 

By  incorporating  the  symmetry  of  the  stress  and  strain  tensors, 
the  quadratic  dependence  of  the  Helmholtz  free  energy  on  strain,  and  the 


Voigt  single  index  notation,  [2]  the  constitutive  equations  reduce  to 
(see  appendix  A) 


(1 


where 


CT'  =  CT  +  M1?.  an 
L  L  .  limn  mn 

1J  iJ 


(2 


In  equations  (1)  and  (2)  the  subscript  L  represents  quantities  which  are 
locally  averaged  at  the  lamina  level.  The  subscripts  i  and  j  range  from 
one  to  six,  the  subscripts  m  and  n  ranges  from  1  to  3i  and  the 
superscript  n  ranges  from  1  to  N,  the  number  of  damage  modes. 

At  this  point  no  further  reductions  can  be  made  to  the  number  of 
unknown  constants  in  equation  (2)  without  specifying  the  specific 
damage  modes  and  material  symmetry. 


Application  to  Matrix  Cracking 

As  discussed  in  the  introduction,  the  potential  of  the 
constitutive  model  will  be  demonstrated  by  considering  the  case  of 
matrix  cracking  in  continuous  fiber  composites.  An  example  of  this 
damage  state  is  shown  in  Fig.  1.  In  order  to  construct  the  proposed 
constitutive  model  for  this  system  we  first  examine  the  response  of  a 
single  ply  subjected  to  transverse  matrix  cracking  as  shown  in  Fig.  2. 
Assuming  that  the  crack  geometry  is  symmetric  about  normals  to  each  of 
the  ply  coordinates,  the  internal  state  variable  associated  with  matrix 
cracking  in  local  ply  coordinates  is  represented  by 


- 
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(3) 


This  implies  that  the  damage  in  a  single  ply  is  defined  by  the  direction 

of  vectors  nc  and  uc  and  the  magnitude  |  a^|  (see  Fig.  2).  It  will  be 

shown  in  a  later  section  that  the  magnitude  of  is  related  to  the 

ij 

surface  area  of  matrix  cracks  per  unit  volume  in  a  ply. 

Note  that  a  second  order  tensor  representation  of  the  internal 
state  variable  is  insufficient  if  the  crack  orientation  within  the  ply 
is  time  dependent.  In  this  case  a  higher  order  tensor  is  required. 
However,  since  the  crack  is  matrix  dominated  and  constrained  by  fibers, 
rotation  is  assumed  to  be  negligible  and  the  second  order  tensorial 
representation  is  considered  adequate  in  this  example. 

For  the  single  damage  mode  of  matrix  cracking  described  in  Fig.  1, 
equation  (2)  reduces  to 


C'  =  CT  +  m!.  a1 
L  .  .  L .  .  ij  mm  tan 

ij  ij 


(4) 


where  a/  represents  matrix  cracking  and  the  above  properties  are 
^mn 

constructed  in  ply  coordinates.  Since  the  only  non-zero  component  of 

a1  is  a.1  ;  equation  (4)  reduces  to 

b  ran  b  2  2 


c;  =  cT  +  a 

L . .  L . .  li 22 

ij  IJ 

where  rePresent  the  non-zero  components  of 

matrix  cracking.  At  this  point  has  Peen 

coefficients  in  equation  (5). 


(5) 

M^.  . .  .  associated  with 
ljkl 

reduced  to  21  unknown 
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Material  symmetries  may  now  be  utilized  to  further  simplify  the 

constitutive  equations.  The  material  in  question  is  assumed  to  be 

initially  transversely  isotropic  in  the  undamaged  state  on  the  local 

scale,  where  the  plane  of  isotropy  is  the  -  x^  plane  shown  in  figure 

2.  In  the  undamaged  state  the  effective  modulus  tensor  c!  is 

Lijkl 

equivalent  to  the  elastic  transversely  isotropic  ply  properties  given  by 


It  can  be  shown  using  irreducible  integrity  bases  [4,5,8]  that 

with  the  inclusion  of  the  damage  tensor  for  tr  ansverse  cracking  the 

effective  modulus  tensor  becomes  orthotropic.  Therefore,  the  damage 

tensor  M,1  is  an  orthotropic  tensor  containing  9  unknown  constants. 
Lij  22 

Thus,  the  complete  constitutive  equation  (5)  (assuming  the  damage  growth 

law  a.  is  known)  requires  the  determination  of  independent  material 
L22 

constants  for  the  undamaged  modulus  tensor  C  ,  and  9  independent 

L.  . 

i  1J 

constants  for  the  damage  tensor,  M 

Lij22 

For  relatively  thin  laminates  it  is  useful  to  apply  the  conditions 
of  generalized  plane  stress  where  the  out-of-plane  shear  stresses  and 
are  neglected.  Applying  these  conditions  to  equation  (1)  and  using 


matrix  notation  results  in 


CL  CT  CT  0 

11  i2  L12 

CL  CL  CL  0 

12  22  23 

c  C  C  0 

12  L23  L22 


M  ^  M*  0 

1122  1222  1322 
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l  1222  2222  1  2322 
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n1322  m2322  n3322 
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where  for  generalized  plane  stress  conditions  there  are  12  unknown 
material  constants  in  equation  (7). 


Laminate  Equations 

To  utilize  single  lamina  equations  to  characterize  the  response  of 
multilayered  laminates,  it  is  necessary  to  globally  average  the  local 
ply  constitutive  equations.  This  is  accomplished  herein  by  imposing  the 
Kirchhoff  hypothesis  for  thin  plates.  This  procedure  is  very  similar  to 
the  method  used  in  Jones  [3].  However,  generalized  plane  strain 
conditions  are  imposed  rather  than  plane  strain  because  this  is 
consistent  with  the  stress  state  in  equation  (7)  (A  detailed  description 
of  the  global  averaging  is  given  in  appendix  B) .  The  resulting 
equations  are  as  follows: 


{N}  =  [A]  {e0} 


l c° }  =  [A  ]  {N  } 


where 


(CL.A  (tk> 

i.l 


%  1  .. 
t'-tvO 


L.-.- 
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and  {.:Q}  contains  the  laminate  midplane  strains  where  k  specifies  the 
ply  and  t  is  the  ply  thickness. 


DEPENDENCE  OF  THE  FREE  ENERGY  ON  CRACK  SURFACE  AREA 


It  was  previously  shown  in  Part  1  [1]  that  the  free  energy,  was 


first  order  in  the  internal  state  variables, 
general 


'ij 


Therefore,  in 


UL"  fij  (EL  ’  V  •  °£.. 
J  kl  13 


(11) 


n 

where  f  is  a  tensor  associated  with  each  specific  mode  of  damage,  n. 

Therefore,  for  the  "fixed  grip"  condition  (e  ■  constants)  and 

Lij 

constant  temperature  conditions,  the  free  energy  is  a  linear  function  of 
the  internal  state  variables. 

The  free  energy,  uc,  depends  on  the  energy  released  during 
cracking  as  well  as  the  total  surface  area  of  cracks  in  the  local 
volume.  For  an  elastic  material  undergoing  stable  self-similar  crack 
growth,  the  available  energy  is  related  to  the  energy  required  for  crack 
extension  by  [6,7] 


U 


c 

L 


PLVL 


CL  S2 


(12) 


where  is  the  local  volume  averaged  strain  energy  release  rate  during 
crack  extension,  and  S?  is  the  total  surface  area  of  matrix  cracks. 
Also,  it  has  been  previously  shown  [1]  that  the  local  energy  is  related 


to  the  total  energy  by 
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The  total  energy  released  during  the  fracture  process  is  obtained  by 
substituting  equation  (13)  into  equation  (12)  and  integrating  in  time. 
This  results  in 


UL(t> 


s9  dt' 


1 _  i 

_*  to  * 

kAVW 


where  t  is  the  time  of  interest,  t'  is  a  dummy  variable  of  integration, 


L,  T  i 


and  S2  is  a  time  dependent  geometric  quantity.  It  should  be  noted  that 
G  may  be  time  dependent  if,  for  example,  cracks  propagate  through 

L« 

heterogeneous  zones  such  as  resin  rich  regions  during  the  damage 


process . 


Since  the  body  is  assumed  to  be  elastic  at  a  constant  damage 


state,  the  strain  energy  release  rate,  G^,  may  be  assumed  to  be  rate 

independent.  Therefore,  equation  (14)  becomes 
S,(t) 


u£(t)  = J  GL(S2)ds  =  g(S2) 


S2(o) 

where  it  is  apparent  that  the  total  energy  is  a  function  of  surface  area 


>s-;v: 


only. 


It  therefore  follows  from  equations  (11),  (13)  and  (15)  that  for 


the  case  of  transverse  matrix  cracking 


ul  =  -L- 

°LVL 

Thus,  finally 


R(S,>  =  ,TL),aL  ' 

12.  LR1  L  L22 


i  =  s(s2)  _  .  K,(s  ) 

22  *)cVc  f  2Z 
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The  precise  nature  of  the  function  g’  is  determined  either  analytically 
(using  a  finite  element  solution)  or  experimentally. 


From  equations  (17)  and  (15),  it  is  seen  that  the  internal  state 
variable  depends  on  the  strain  energy  release  rate,  GL-  For  the  "fixed 
grip"  condition,  the  strain  energy  release  rate  is  given  by  [7] 

G  =-i(^T)  ,  U8) 

B  da 

where  B  is  the  crack  depth,  and  UT  is  the  elastic  strain  energy.  For  a 
uniaxial  stress  state,  the  elastic  strain  energy  is  defined  as 


ut  E/[/and€n]  dv  •  <19> 

where  °and  are  the  uniaxial  stress  and  strain  states, 

respectively,  and  V  is  the  volume  of  interest.  For  a  linear  elastic 
material,  the  term  in  the  brackets  of  equation  (19)  can  be  integrated 
and  equation  (19)  becomes 

UT  ■/ [|  E11  Ell]  dV  <20) 

where  E1 1  is  the  laminate  stiffness  (A^/t)  in  the  x-axis  direction. 
For  the  uniaxial  fixed  grip  condition,  e  does  not  vary  over  the 
volume.  Therefore,  equation  (20)  becomes 


UT  "I 
2 

v  En  En 

• 

(21) 

Finally, 

substituting 

equation 

(21)  into  equation  (18) 

and 

recalling 

that 

is  constant 

for  the 

"fixed  grip"  condition, 

the 

following 

expression  for  strain 

energy  release  rate  is  obtained: 

gl  “  - 

1  Vc2  ')E11 
2B  11  3a 

(22) 

viviviv 


12 


i 


This  result  implies  that  only  the  laminate  modulus,  E^  changes  with  the 
formation  of  new  crack  surfaces. 

The  above  result  may  be  considered  intuitively  obvious  because  E j 
is  defined  as  the  load  in  the  x^-direction  divided  by  the  associated 
displacement  in  the  x  j-direction.  Since  the  displacement  is  held 
constant,  the  load  will  decrease  as  crack  extension  occurs  and  the 
structure  becomes  more  compliant.  This  can  be  shown  to  be  analogous  to 
[7] 

Gr  =  li  9£  *  (23) 

2B  3A 

where  C  is  the  specimen  compliance  (1/E^).  . 

As  an  example,  consider  a  [0/90]^  laminate  with  matrix  cracking  in 
the  90°  plies  only.  The  average  stiffness  is  given  by 

E11  ■  <E0+  E90>  n  •  <24) 

where  EQ  and  EgQ  are  the  moduli  in  the  global  x1  direction  (0°  fiber 
orientation)  of  the  0°  and  90°  plies,  respectively.  Assuming  that  no 
cracks  form  in  the  outer  0°  plies,  equation  (22)  becomes 

G  =  -  J_  Ve2  3E90  .  (25) 

L  4B  3a 

This  result  implies  that  only  changes  in  the  modulus  of  the  constrained 
90°  plies,  where  matrix  cracks  are  assumed  to  occur,  effect  the  value  of 
the  strain  energy  release  rate  for  increasing  crack  formation. 


Determination  of  Damage  Modulus  Tensor 

The  material  constraints  used  in  M 


1 

■U  22 


can  be  determined  using 


equation  (5) 


provided  that  C,  ,  C.'  ,  and  a}  are  known.  In  reality  it  is  Ct* 

Lij  Lij  L22 

that  one  seeks  to  predict;  however,  it  is  possible  to  analytically  or 

'  1 

experimentally  determine  C  for  prescribed  values  of  a.  .  The  terms 

Lij  L22 

in  C  are  the  orthotropic  material  properties  determined 

Li  •  • 

13 

experimentally  from  tests  conducted  on  undamaged  composite  lamina. 

f 

In  order  to  determine  the  terms  in  C.  ,  it  is  necessary  to 

Li  ,  . 

determine  the  effective  orthotropic  properties  of  a  single  ply  for  a 

given  value  of  a/  .  The  effective  orthotropic  properties  due  to 
L22 

matrix  cracking  are  determined  by  modeling  the  cracky  as  internal 

boundaries  and  solving  the  resultant  boundary  value  problem  using  a 

» 

finite  element  continuum  model.  Alternatively,  C  may  be  determined 

Lij 

experimentally  where  methods  exist  to  determine  appropriate  stiffness 
components . 


RESULTS 

This  section  presents  the  results  obtained  from  the  finite  element 

continuum  model  which  was  used  to  predict  the  damage  modulus  tensor  and 

the  globally  averaged  laminate  equations  relating  stress  and  strain. 

These  results  demonstrate  the  ability  of  the  model  to  predict  all 

components  of  the  effective  stiffness  due  to  matrix  cracking  for  a 

[0,90  ]  laminate.  The  experimental  data  used  in  this  analysis,  shown 
j  o 

in  Table  1,  were  obtained  from  reference  9. 

The  first  series  of  tests  were  conducted  to  determine  the  strain 
energy  release  rate  and  the  damage  modulus  tensor  m' ?  for  [0,90]  and 

1  J  3 

[90,90]s  laminates  using  the  finite  element  continuum  model  with  various 


crack  densities  in  the  center  90°  ply.  The  finite  element  me3h  used  for 
the  analysis  is  shown  in  Fig.  3.  The  effect  of  different  constraining 
layers  on  the  strain  energy  release  rates  calculated  from  the  finite 
element  results  are  shown  in  Fig-  4.  Also  shown  in  Fig.  4  is  the 
theoretical  result  for  repeated  cracks  in  an  isotropic  infinite 
domain  [10].  This  result  is  intended  only  to  verify  the  cubic  nature  of 
the  strain  energy  release  rate  calculated  from  the  finite  element  model. 
For  the  current  state  of  model  development,  only  the  [0,90]a  relation 
for  the  energy  release  rate  was  used.  It  should  be  noted  that,  in 
calculating  the  energy  release  rate,  the  derivative  of  the  change  in 

l 

stiffness  with  respect  to  crack  length  was  approximated  with  a  first 
order  finite  difference  relation. 

The  next  step  was  to  determine  an  approximate  relation  between  the 

strain  energy  release  rate  and  the  total  crack  length.  A  least  squares 

cubic  fit  of  the  curve  for  the  [0,90]  laminate,  shown  in  Fig.  4  was 

s 

obtained  and  is  given  as  follows: 

G(a)  =  3.88690  -  0.280162  a  -44.11385  a2  +  69.22121  a3,  (27) 

where  the  nonlinear  terms  are  believed  to  be  caused  by  boundary  effects 
and  crack  interaction.  Substituting  equation  (27)  into  equation  (15) 
results  in  an  expression  for  the  internal  state  variable  as  follows: 

aL  =  3.8869a  -  0.140081  a2  -  14.70462  a3  +  17.305303  a^  .  (28) 

Tne  values  for  the  damage  modulus  tensor  can  now  8e 

determined  from  the  finite  element  results  for  the  [0,90]^  laminate.  In 
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...... 


order  to  evaluate  the  constants  in  M.j ^ ,  a  particular  crack  density  is 

selected  and  the  finite  element  results  for  this  case  are  utilized  in 

equation  (25).  These  results  are  then  used  as  follows:  1)  compare 

damage  model  predictions  for  the  axial  stiffness  loss  in  the  [0,90]s 

laminate  to  the  finite  element  results,  and  2)  compare  the  damage  model 

predictions  for  the  axial  stiffness  loss  in  a  [0,90  ]  laminate  to  both 

5  o 


finite  element  predictions  and  experimental  results. 

The  results  for  case  one  are  shown  in  Fig.  5.  It  can  be  seen 
that  the  damage  model  underpredicts  the  stiffness  loss.  This  result  is 
expected  since  laminate  analysis  does  not  account  for  interlaminar 
shearing  which  would  effectively  reduce  the  axial  stiffness. 
Furthermore,  the  strain  energy  release  rate  predicted  by  the  finite 
element  continuum  model  is  conservative  due  to  the  inherent  "over 
stiffness"  problems  associated  with  the  finite  element  model,  especially 
since  constant  strain  triangular  elements  were  incorporated.  This  same 
argument  will  hold  true  when  the  finite  element  model  results  are 
compared  to  experimental  data. 

The  normalized  axial  stiffness  results  for  case  2  ( [ 0 , 90_ ] 

d  s 

laminate)  are  shown  in  Fig.  6.  As  discussed  previously,  the  results 
predicted  from  the  damage  model  as  well  as  finite  elements  underestimate 
the  stiffness  loss  in  the  laminate.  Nevertheless  the  good  comparison  of 
the  damage  model  to  the  experimental  results  can  only  be  encouraging 
since  the  only  experimental  data  used  were  the  undamaged  stiffness 
properties.  The  inclusion  of  experimental  data  for  the  damaged 
stiffness  properties  and  a  more  refined  energy  release  rate  prediction 
can  only  help  decrease  the  current  differences  between  model  predictions 
and  experimental  results.  A  phenomenological  explanation  for  the 
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6.  Normalized  Axial  Stiffness  Predictions  for 
a  [0,90j]s  E-Glass  Epoxy  Laminate. 
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underpredi cted  stiffness  loss  can  be  seen.  If  one  examines  a  schematic 

diagram  of  a  typical  saturation  crack  pattern  for  the  [0,90  ]  laminate 

j  s 

[9,11]  as  shown  in  Fig.  7a.  It  is  observed  that  a  network  of  branched 

or  curved  cracks  exist  in  addition  to  a  more  uniform  spacing  of 

transverse  cracks.  This  suggests  that  the  surface  area  of  transverse 

matrix  cracks  may  actually  be  greater  than  that  assumed  by  the  damage 

model  as  shown  in  Fig.  7b.  Since  the  internal  state  variable  depends 

directly  on  the  surface  area  of  cracks,  one  concludes  that  by  allowing 

for  the  increased  surface  area  due  to  branching  or  curvature  then  the 

model  results  would  be  in  better  agreement  with  the  experimental  data. 

This  suggests  that  it  may  be  preferrable  to  obtain*  all  material 

parameters  from  experimental  tests  rather  than  finite  element  results. 

In  Fi;.; .  8  the  remaining  components  of  the  effective  modulus 

tensor  are  predicted  using  the  laminate  equations.  It  can  be  seen  from 

the  figure  that  the  moduli  in  the  y  and  z  directions  do  not  vary 

significantly  for  transverse  matrix  cracking.  However,  it  can  be  seen 

that  the  Poisson's  ratios  v  and  v  ,  as  well  as  the  shear  modulus 

xy  xz 

G  ,  vary  significantly.  Also  shown  in  Fig.  8  are  results  obtained 

xy 

from  the  finite  element  continuum  model  for  The  good  agreement 

demonstrated  here  supports  the  ability  of  the  model  to  predict  all 
components  of  the  effective  modulus  tensor. 

CONCLUSION 

In  part  I  constitutive  relations  were  developed  to  account  for 
various  forms  of  damage  in  continuous  fiber  composites  by  utilizing  the 
concept  of  internal  state  variables.  These  constitutive  equations  have 
been  simplified  to  a  usable  form  by  imposing  symmetry  constraints .  The 
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Fig.  8.  Effective  Orthotropic  Property  Predictions  for 
[0,903Js  E-Gla  ss  Epoxy  Laminate. 
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potential  of  the  resulting  constitutive  model  was  demonstrated  by 
considering  the  case  of  matrix  cracking  only.  Results  were  obtained  for 
[0,90]  and  [0,90,]  laminates  with  matrix  cracking  in  the  90°  plies. 
Although  some  of  the  material  constants  used  to  characterize  the  local 
ply  constitutive  equations  were  determined  analytically,  model  results 
compare  favorable  to  experimental  and  finite  element  analyses  for 
specified  values  of  crack  densities.  It  is  anticipated  that  when  the 
material  constants  used  in  the  model  are  determined  experimentally,  the 
results  should  be  notably  better  due  to  the  inability  of  the  finite 
element  results  to  account  for  unusual  crack  geometry.  Finally,  the 
model  was  shown  to  be  capable  of  predicting  all  components  of  the 
damage  induced  effective  modulus  tensor.  However,  the  accuracy  of  this 
capability  remains  to  be  verified  through  experimental  correlation. 

The  objective  of  this  research  has  been  threefold:  1)  to  develop 
a  model  for  predicting  the  constitutive  behavior  of  a  laminate  given  the 
mode  and  extent  of  damage,  2)  to  develop  growth  laws  for  each  mode  and 
the  extent  of  damage,  and  3)  to  incorporate  both  1)  and  2)  into  a  usable 
form.  Although  the  model  has  been  demonstrated  for  given  states  of 
matrix  cracking,  parts  2)  and  3)  require  further  research  and  are  indeed 
the  most  difficult  link  in  the  model  development  since  they  will  require 
extensive  correlation  with  experiment  as  well  as  finite  element 
analysis . 
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Appendix  A:  Application  of  Symmetry  Constraints 


The  damage-dependent  constitutive  model  (equations  (61*)  through 
(67)  of  Part  I  [1])  is  defined  as  follows: 


R  *  T 

a  =  a  +  C  (e  -  e  ) 

L .  .  L  L .  ..  ,  L,  ,  kl 

ij  ij  ijkl  kl 


where  o.  -  local  stress  tensor 
Lij 


local  strain  tensor 


o,  S  0,  (B,  +  I?...  ) 

llj  L  Lij  ljkl  lkl 


Residual  stress  in  the  absence  of  strain  and  temperature 
change 


Ct  2  pr  <Ct  +  ) 

L..,,  L  L.,.,  ljklmn  L 

ijkl  i.]kl  J  mn 


effective  modulus  tensor 


C  =  undamaged  modulus  tensor 
ijkl 

4  =  4'  °L  <EL..  -4kl“U  “l> 

ijkl  ij  J 


thermal  strain  tensor 


internal  state  variable  for  damage  modes 


change  in  temperature 


atl  e 


and  B  ,  E  &  ... 

L„  ijkl  ijklmn,  ijkl 


are 


tensorial  material  constants  as  required  by  the  Taylor  Series  expansion. 

For  demonstration  purposes,  the  residual  stress  terms  are 
neglected  and  isothermal  behavior  is  imposed  on  equations  (la)  through 
(4a),  yielding 


°L  ~  er 

ij  ijkl  kl 


where  the  effective  modulus  tensor  is  given  by 


(5a) 


n  \ 

aT  ) . 


p.  (C,  +  M  '  .  „T 

L  L .  ..  .  ijklmn  L 

ijkl  mn 


(6a) 


Note  that  M. .  is  a  sixth  order  tensor  with  729  coefficients  for  each 
ljklmn 

value  of  n.  It  is  assumed  here  that  the  constitutive  equations  given  by 
(5a)  are  statistically  homogeneous.  Therefore,  the  conditions  of  stress 
and  strain  symmetry  can  be  applied  to  equations  (5a)  and  (6a)  to  obtain 


CL  =  CL 
ijkl 


jikl 


CL  =  CL 
ijkl  ijlk 


(7a) 


Using  equations  (6a)  and  (7a)  it  can  be  shown  that 


M*1  =  Mn 

ijklmn  jiklmn’  ijklmn  ijlkmn 


(8a) 


Using  the  Helmholtz  free  energy  and  the  quadratic  dependence  of  the 


effective  modulus  tensor  on  strain,  it  can  be  shown  that 


Substitution  of  equation  (6a)  into  (9a)  will  result  in 


MV ,  =  M.r',  .  .  (10a) 
ljklmn  klijmn  v  a' 

The  above  symmetry  constraints  reduce  the  number  of  coefficients  to  189 
terms  in  equation  (10a)  for  each  value  of  n.  It  is  most  convenient  at 
this  point  to  reindex  the  constitutive  tensors  using  the  Voigt  notation 
[2]  where 


o 

11 

°11 

°4  °23  “  °32 

II 

<M 

o 

°22 

°5  =  °13  =  o31 

°3  = 

°33 

°6  =  °12  °21 

and 


e1  C1 1 

e4 

=  2 

e23  =  2  e32 

£_  =  £*- 

£r 

=  2 

C  .  «=  2 

(12a) 

2  22 

5 

13  31 

G3  =  e33 

e6 

=  2 

p  s  2  £ 

e1 2  21 

Using  the  contracted  notation,  equations  (5a)  and  (6a)  can  be  written  as 


a  =  C  e. 
L .  L »  .  3 

i  ij 


where  i  and  j  range  from  1  to  6  and 


(13a) 


C.  =  C  +  M,.  an 
L,  .  L.  .  ljmn  mn 
ij  ij 


(14a) 

Note  that  m  and  n  range  from  1  to  3  and  n  ranges  from  1  to  N,  where  N  is 


the  number  of  damage  modes. 


Appendix  B 


The  values  of  generalized  plane  strain  are  given  by 
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(1b) 


where  the  superscript  "o"  denotes  the  midsurface  strains  and  the  < 

matrix  denotes  the  midsurface  curvatures.  Under  the  condition  of 

generalized  plane  strain  there  is  no  warping  allowed  out-of-plane,  which 

implies,  that  k  »  0.  Thus  the  generalized  plane  strain  equation 
z 

becomes 


c 

c°  ' 

t* 

X 

t. 

X 

e 

y 

z 

>  =  < 

z° 

y  > 

e°  ( 

z 

z 

z 

xy 

£° 

xy 

+ 


t 


(2b) 


It  is  now  assumed  that  no  moments  or  curvatures  are  imposed  and 
that  all  laminates  studied  are  symmetric.  Therefore,  in  order  to 
determine  the  resultant  forces,  it  is  necessary  only  to  integrate  the 
given  stress  state  as  specified  in  equation  (7)  over  the  laminate 


thickness  to  obtain 


where  t  is  the  total  thickness  of  the  laminate. 

Substituting  equations  (13a)  and  (2b)  into  (3b)  for  the  case  where 
there  are  no  rotations  results  in 


{N}  = 


J  CC 

-t/2 


t/2 

CC‘] 


{e°}  dz 


(4b) 


where  {N}  denotes  the  force  resultants,  [C'l  is  the  transformed 

Li 

effective  stiffness  matrix,  and  {e°}  represent  the  mid-surface  strains. 
Note  that  since  transverse  cracks  go  completely  through  the  thickness  of 
the  cracked  plies  the  stiffness  is  assumed  to  be  spacially  constant 
through  the  thickness.  Therefore,  equation  (4b)  can  be  written  as 

{N}  =  Z  [C^]k  (ZR  -  Zk-1)  {c°}  ,  (5b) 

k=l 


where  k  specifies  the  ply  and  Zk  -  Zk_^  is  the  thickness  of  each  ply. 


One  can  define 


where  A  represents  the  globally  averaged  stiffness  matrix.  Thus  the 
globally  averaged  constitutive  equations  become 

{Nj  =  [A]  {e0}.  (7b) 

Experimental  testing  is  often  conducted  on  uniaxial  testing 
machines  in  which  the  applied  force  resultants  are  input  and  the  strains 
are  experimentally  determined  output.  Therefore,  at  times,  it  is  more 
convenient  to  express  the  strains  in  terms  of  the  applied  force 


resultants  as  follows. 


